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FOREWORD 


This  report  discusses  a  three-degree-of-freedom  dynamic 
wind  tunnel  test  program  using  slender  cones  with  and  without 
nose- tip  ablation.  These  tests  were  carried  out  In  the  NSWC 
Hypersonic  Tunnel  using  a  gas-bearing  support  designed  and 
fabricated  specifically  for  this  program.  A  data  reduction 
program,  also  formulated  for  these  tests,  provided  static  and 
dynamic  stability  derivatives  from  measurements  of  bl-planar 
displacement  angles.  Comparisons  are  made  where  possible  with 
results  from  previously  conducted  static  tests.  A  significant 
result  of  these  tests  was  a  demonstration  of  the  dynamic  un- 
stabllltles  associated  with  the  presence  of  nose-tip  ablation. 
Also,  It  was  shown  tiiat  during  ablation  there  Is  continual  and 
significant  forward  motion  of  the  static  center  of  pressure. 
This  test  program  was  sponsored  by  the  Reentry  Aerodynamics 
program. 
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1.0  INTROPOQTlOH 

In  conventional  wind  tunnel  testing  the  model  is  attached  to  a  rigid  struc- 
v  ture,  often  a  cantilever  beam.  This  support  is  called  a  "wind  tunnel  balance"  since 

an  integral  part  of  the  structure  is  a  load  measuring  unit.  A  measurement  of  die 
forces  of  constraint  is  assumed  equivalent  to  a  measurement  of  the  aerodynamic 
loads.  By  rotating  the  model  and/or  support  structure  relative  to  the  flow  vector, 
the  dependency  of  load  on  angle,  can  also  be  measured.  While  there  is  much  to 
recommend  testing  of  rigidly  constrained  models,  such  testing  cannot  measure  the 
aerodynamic  effects  caused  by  the  rate  of  change  of  angular  variables.  The  three- 
degree-of-freedom  gas  bearing  support  was  designed  to  permit  the  measurement  of 
time  dependent  aerodynamic  loads.  The  model  is  supported  by  a  cantilever  beam  with 
the  attachment  point  at  the  model's  center  of  gravity.  The  support  permits  re¬ 
stricted  angular  motion  in  pitch  and  yaw  {about  7.5  degrees)  and  unrestricted 
angular  motion  in  roll. 

While  physical  realism  may  be  enhanced  by  the  testing  of  unconstrained  models, 
the  determination  of  the  aerodynamic  loads  is  no  longer  a  primary  measurement  as  it 
is  in  constrained  model  testing.  The  primary  measurement  is  model  angular  dis¬ 
placement  from  a  null  or  trim  condition.  In  the  three-degree-of-freedom  gas 
bearing  support  these  angles  are  roll,  pitch  (angle  of  attack),  and  yaw  (angle 
of  sideslip).  The  task  facing  the  analyst,  then,  is  to  deduce  what  die  loads 
must  have  been  during  the  test  in  order  to  produce  die  observed  motion.  The  pro¬ 
cedure  then  involves  postulating  or  modeling  the  flow-body  interaction.  This  is 
equivalent  to  postulating  the  form  of  relevant  differential  equation (s). 

In  the  present  application  the  differential  equation  describes  the  external 
flow  of  a  compressible  fluid  about  a  n on-permeable  structure.  The  fluid  is 
*  regarded  as  being  everywhere  continuous  and  the  component  of  the  fluid  velocity 

normal  to  the  surface  of  the  structure  is  everywhere  zero.  If  the  differential 
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•qttttloB  Is  linear  then  the  coefficients  may  be  regarded  as  "constants  of  propor¬ 
tionality"  between  angular  displacement  (and  its  derivative)  and  the  applied 
moment,  the  recordings  of  angular  displacement  (primary  measurements)  are  in  a 
sense  "solutions"  of  the  differential  equations;  the  goal  of  data  reduction,  then, 
is  to  determine  the  numerical  values  of  the  coefficients  that  provide  (according  to 
some  criterion)  the  "best  fit"  to  the  observed  (measured)  angular  motion.  The 
criterion  that  will  be  used  here  for  "best  fit"  will  minimise  the  square  of  the 
"error"  or  difference  between  the  assumed  function  and  the  data  points. 

The  emphasis  of  this  report  will  be  on  the  use  of  the  Three-Degree-of-Freedom 
Gas  Bearing  (TDGB)  rather  than  on  its  mechanical  design.  However,  a  description 
of  the  — chemical  features  along  with  annotated  assembly  drawings  and  photographs 
will  be  Included.  The  readout  or  date-gathering  techniques  associated  with  this 
TDGB  can  never  be  regarded  as  complete,  since  technological  improvements  will 
suggest  alterations  from  time- to- time. 

By  way  of  supporting  claims  of  utility,  some  wind  tunnel  test  results  will 
also  be  included.  Non-dimensional  aerodynamic  coefficients  will  be  given  for  both 
ablating  and  non-ablating  configurations. 

2.0  BEARING  DBSIGH/Q0N8TKPCTI0N 

The  Three-Degree-of-Freedom  Gas  Bearing  (TDGB)  permits  unlimited  model 
freedom  in  roll  and  restricted  freedom  in  angles  of  attack  and  sideslip.  It  is 
certainly  possible  to  provide  the  same  degrees  of  freedom  with  a  global  system 
using  mechanical  bearlnga.  However,  In  the  preaent  application  where  heat  loads 
are  modest,  a  gas  bearing  is  particularly  attractive  because  of  its  low  and  con¬ 
sistent  friction.  However,  the  gas  bearing  does  require  a  continuous  supply  of 
gas  (air  in  prssant  case)  and  as  it  turns  out,  some  control  of  the  temperature  of 
the  gas  entering  the  region  of  the  bearing  sphere  is  necessary  to  maintain  sphere/ 
'socket  clearances.  The  bearing  air  also  must  exit  the  model's  base.  However  die 
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mass  flow  is  so  small  that  exhaust  air  has  no  significant  effect  on  model  motion. 

The  sphere/ cone  model  is  shown  attached  to  the  gas  bearing  and  mounted  in  the 
NSWC  Hypersonic  Wind  Tunnel  in  Figure  1.  This  model  is  a  7°  cone  (half-angle) 
with  at  22  percent  spherical  bluntness.  All  aerodynamic  data  presented  in  this 
report  were  ootained  using  this  model  (see  section  (4)). 

The  gas  bearing  assembly  and  angular  readout  is  shown  in  Figure  2  .  This 
support  and  readout  system  is  intended  for  conical  models.  The  model  mounts 
directly  to  the  assembly  which  in  turn  is  attached  to  the  cantilever  or  "sting" 
support. 

Figure  3  is  an  exploded  view  of  the  assembly  shown  in  Figure  2.  The  front 
and  aft  ball  housings  form  the  stator  of  the  2.5  inch  diameter  sphere  also  called 
the  "socket."  The  inner  surface  of  the  housings  is  molded  epoxy  (DEVCON) .  The 
clearance  tolerance  between  the  sphere  and  the  epoxy  socket  is  on  the  order  of 
0.001  inch.  The  method  of  forming  this  socket  and  maintaining  this  tolerance 
is  interesting.  When  the  sphere  is  fabricated  a  second  and  nearly  identical  sphere 
is  also  fabricated.  This  second  sphere  is  then  nickel  plated.  The  epoxy  socket 
is  then  molded  to  this  plated  sphere.  The  unplated  sphere  then  fits  the  socket 
with  the  required  clearance.  Sphere-socket  clearance  may  vary  somewhat  so  the 
plating  thickness  may  vary  a  few  thousandths  of  an  inch.  However,  more  important 
than  the  clearance  itself  is  consistency  in  clearance  all  around  the  Dearing 
surface.  The  only  difficulties  experienced  with  the  bearing  clearance  were  caused 
by  cooled  air  (due  to  expansion)  which  contracted  the  epoxy  after  prolonged  bearing 
use.  Heating  the  air  prior  to  entry  at  the  bearing  surface  totally  eliminated  this 
problem. 

All  components  aft  of  the  bail  housing  are  associated  with  tne  pitch/yaw/roll 
fiber  optic  readout  system.  This  system  was  not  entirely  satisfactory  although 
some  angular  measurements  using  the  fiber  optic  system  were  used  for  da^a 
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FIGURE  1  SPHERP/CONE  MODEL  MOUNTED  IN  THE  NAVAL  SURFACE  WEAPONS  CENTER  HYPER¬ 
SONIC  WIND  TUNNEL. 
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FIGURE  3  GAS  BEARING  ASSEMBLY/READOUT  COMPONENTS. 


12 


NSWC  TR  81-491 


reduction.  Roll  or  spin  rates  measured  using  this  system  provided  a  check  on  the 
Independent  measurements  that  led  to  computations  of  the  two  oscillatory  fre¬ 
quencies  (see  equation  (11c)). 

The  fiber  optic  system  readout  will  not  be  discussed  in  great  detail  since  it 
is  to  be  subjected  to  modification  or  replacement.  In  principle  the  fiber  optic 
readout  is  fairly  straightforward,  although  in  practice  fiber  breakage  limited  the 
angular  coverage. 

The  fiber  optic  light  guides  are  0.020  inch  diameter.  These  guides  are 
located  at  discrete  angular  positions.  A  light  source  outside  of  the  tunnel  is 
conducted  by  the  guides  through  the  sting  and  inside  the  bearing.  The  angles  of 
attack  (pitch)  and  sideslip  (yaw)  are  read  approximately  every  half  degree;  the 
roll  angle  is  read  every  22.5  degrees.  As  the  model  rotates  a  mask  covers  and 
uncovers  the  guides  from  the  external  light  source.  As  a  guide  is  uncovered  the 
light  crosses  a  0.050  inch  gap  to  enter  another  guide  directly  opposite  the  light** 
carrying-guide.  The  mask  sheath  (attached  to  the  model  to  share  its  angular 
rotation)  moves  within  this  gap. 

The  light  from  the  pickup  guide  then  goes  to  a  photocell  where  it  is  converted 
to  an  electrical  signal.  A  total  of  five  channels  are  required,  two  for  pitch, 
two  for  yaw  and  one  for  roll.  After  these  signals  have  been  recorded  (analog)  on 
magnetic  tape,  they  are  digitized.  The  digitized  signal,  together  with  a  "readout 
algorithm"  and  calibration  provides  the  angular  record,  l.e.  degrees  (In  yaw  and 
pitch)  versus  time.  In  the  next  section  the  discussion  will  cover  how  this 
angular  record  be  used  to  calculate  the  characteristic  frequencies  and  damping 
exponents  end  ultimately  the  non-dimensional  aerodynamic  derivatives. 

Figure  4  provides  a  schematic  of  the  more  salient  features  of  the  bearing. 

It  will  be  noted  that  the  sting  is  not  solid  but  contains  various  passages  for  the 
light  guides  as  well  as  conduits  for  bearing  and  control  air.  There  are  five 
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separate  air  supply  passages:  (a)  bearing,  (b)  kicking  jet,  (c)  centering 
fingers- forward,  (d)  centering  fingers  retract,  (e)  turbine  (see  Figure  5). 

Bearing  air  Is  admitted  into  the  sting  at  about  375  to  400  psi.  As  pointed 
out  earlier,  bearing  air  should  be  heated  so  that  there  is  no  significant  contrac¬ 
tion  of  the  epoxy  socket.  Setting  an  exact  temperature  is  not  feasible,  nor  Is  it 
necessary:  the  air  after  expansion  In  the  bearing  must  be  equal  to  or  somewhat 
greater  than  the  socket  temperature.  The  bearing  air  supply  enters  the  bearing 
sphere  as  Indicated  and  exhausts  through  eight  holes — four  forward  and  four  aft. 
This  air  Is  then  carried  by  the  bearing  exhaust  lines  In  the  sphere  to  the  rear  of 
the  bearing  where  It  enters  the  model  Interior.  Bearing  exhaust  air  passes  near 
the  fiber  optic  readout  system,  partially  contributing  to  the  cooling  requirements. 
It  will  be  noted  that  there  is  also  water  cooling  provided  to  the  readout  system. 
The  air  finally  exits  the  base  of  the  model. 

The  second  air  source  is  the  air  jet.  The  air  jet  is  simply  a  high  pressure 
jet  of  air  (about  850  psi)  which  when  applied  to  one  side  of  die  model  Impulsively 
gives  the  model  an  Initial  disturbance  in  angle  and  angular  rate.  Because  the 
model  is  spinning,  a  mechanical  initiator  would  be  unduly  complicated.  The  air 
jet  has  been  found  to  be  an  entirely  reliable  and  effective  means  for  initiating 
model  motion. 

The  third  and  fourth  air  supplies  move  the  piston  and  the  centering  fingers 
forward  and  aft.  Air  at  about  50  psi  enters  the  base  through  the  conduit  indicated 
in  Figure  5a  (and  the  enlargement  In  ELgure  5b).  The  piston  moves  forward, 
forcing  the  centering  fingers  into  the  O-ring.  The  turbine  then  imparts  spin  to 
the  model  while  at  the  same  time  holding  the  model  at  null  In  both  orthogonal 
aerodynamic  angles.  When  the  desired  spin  rate  is  achieved,  air  is  simultaneously 
bled  from  the  region  aft  of  the  piston  and  applied  to  the  forward  piston  face  by 
the  aft  piston  air  supply.  The  piston  and  centering  fingers  are  driven  aft.  The 
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no  del  Is  now  available  for  angle  of  attack  Initiation  by  the  kicking  jet  described 
earlier. 

The  final  air  supply  drives  the  turbine.  Air  is  admitted  to  the  turbine  at 
200  to  250  psi.  The  turbine  does  not  contain  blades  as  night  be  expected  but 
rather  consists  of  18  equally  spaced  holes  or  noscles  drilled  in  the  piston  forning 
a  ring  around  the  piston.  Two  of  these  nozzles  are  shown  in  Figures  4  and  5. 

The  axis  of  these  nozzles  are  directed  such  that  the  effluent  vector  has  a  com¬ 
ponent  tangential  to  the  piston  surface.  When  the  piston  is  pressured  in  a 
positive  direction,  the  centering  fingers  are  forced  against  the  0-ring  attached 
to  the  nodal.  Air  imparted  to  the  turbine  then  spins  the  model  while  at  the  same 
time  holding  the  model's  centerline  directly  into  the  wind  tunnel  flow  vector. 

Once  the  model  has  been  spun  to  the  desired  spin  rate  the  fingers  are  retracted, 
turbine  air  supply  terminated  and  the  kicking  jet  initiated.  At  this  point  the 
nodel  is  assumed  undergoing  spin/yaw/pitch  notion  that  is  satisfactory  for  data 
acquisition. 

Since  the  model  must  be  dynamically  balanced  for  proper  angular  response, 
there  is  provision  for  ensuring  that  the  center  of  mass  is  at  the  point  of  rotation. 
Fore  and  aft  movement  of  the  center  of  mass  is  controlled  by  the  location  of  the 
HEAVIMET  (tungsten  alloy)  balancing  slug  seen  in  Figure  5.  Since  the  nodel  is 
configurationally  a  body  of  revolution,  nost  of  the  components  also  have  an  axis 
of  symmetry  and  are  shaped  on  a  lathe.  However,  internal  screw  holes  and  other 
irregularities  can  place  the  center  of  mass  off  the  configurational  axis  of  sym¬ 
metry.  Consequently  a  balancing  ring  (located  at  the  rear  of  the  model)  can  be 
machined  in  accordance  with  balancing  machine  instructions  to  bring  the  center  of 
mass  laterally  to  the  axis  of  configuration  symmetry. 

Because  the  dynamic  model  is  used  in  a  flow  environment  where  the  stagnation 
temperature  may  be  1700  degrees  Fahrenheit,  prolonged  usage  might  cause  the 
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temperature  In  the  vicinity  of  the  fiber  optic  light  guides  to  exceed  175  degrees 
Fahrenheit.  Thus,  water  cooling  is  provided  by  a  paasage  through  the  center  of  the 
sting.  In  the  present  set  of  tests  water  cooling  was  probably  unnecessary  since 
bearing  effluent  air  provided  sufficient  cooling  for  the  test  durations  (less  than 
60  seconds). 

The  preceding  description  covered  the  essentials  of  the  Model's  Mechanical 
design  features.  The  fiber  optic  system  functioned  only  intermittently  and 
hence  cinematic  coverage  was  used  for  all  data  reduction.  The  only  other  diffi¬ 
culty  experienced  was  binding  (intermittent  stator/rotor  contact)  in  die  gas  bear¬ 
ing  after  prolonged  use.  This  problem  was  alleviated  by  heating  the  bearing  air 
to  about  50  degrees  over  ambient.  Tests  for  maximum  sustainable  loads  on  the  gas 
bearing  have  not  been  made;  however,  it  has  shown  the  capability  of  sustaining  a 
model  weight  of  35  pounds  and  an  axial  aerodynamic  load  of  over  100  pounds. 


3.0  DATA  REDUCTION 


In  this  section  die  relevant  differential  equation  of  no  del  notion  and  the 
solution  of  this  equation  will  be  given.  The  least  squares  date  reduction  tech¬ 
nique  will  be  developed  and  the  coaputer  program,  formulating  this  technique  will 
be  given. 

3.1  DIFFERENTIAL  EQUATION  OF  MOTION.  The  wind  tunnel  nodal  has  been  indi¬ 
cated  to  undergo  angular  notion  about  three  autually  perpendicular  axes  while 
innersed  in  a  continuous  newtonian  nedlun.  The  situation  then  Is  sinllar  to 
classical  aeroballistlc  notion  with  two  restrictions.  In  classical  aeroballls tics 
the  nodel  (or  projectile)  aovea  through  an  essentially  stationary  nedlun  with 
permissible  motion  about  and  along  three  autually  perpendicular  axes.  The 
Gallilean  transformation  from  a  moving  nodel/ 8 tatlonary  nedlun  to  a  linear  con¬ 
strained  model/moving  medium  does  not  change  the  physics  of  aeroballls tics. 
However,  constraining  the  model  against  linear  motion  does  mean  that  only  the 
equations  of  angular  motion  need  be  considered:  the  equations  of  linear  notion 
are  identically  zero. 

There  is  some  controversy  associated  with  the  origins  of  the  equations  of 
aeroballistlc  notion  known  as  "tricyclic  theory."  Reference  (1)  is  often  con¬ 
sidered  the  seminal  paper  in  the  statement  of  the  tricyclic  theory;  reference  (2) 
is  also  a  broad  statement  of  the  theory  and  contains  many  extensions  and  applica¬ 
tions  to  the  ballistic  range.  However,  a  possibly  more  readable  paper  than  either 
of  these  and  one  that  will  be  used  here  as  the  basic  reference  is  reference  (3) 
by  H.  R.  Vaughn. 


^Nicollades ,  J.D.,  On  the  Free  Flight  Motion  of  Missiles  Having  Slight  Con¬ 
figurational  Asymmetries,"  Ballistic  Research  Laboratories  Report  No.  858, 
June  1953. 


Sfurphy,  C.H.,  "Free  Flight  Motion  of  Syjpetrlc  Missiles,"  Ballistic  Research 
Laboratories  Report  No.  1216,  July  1963. 

^Vaughn,  H.R. ,  "A  Detailed  Development  of  the  Tricyclic  Theory,"  Sandie 
Laboratories,  SC-M-67-2933,  February  1968. 
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Relying  on  the  development  contained  In  reference  (3)  it  Is  possible  to  show 
that  the  angular  notion  about  two  orthogonal  and  wind  tunnel-fixed  axes  Is  described 
by  the  equations: 

Io  -  p8I  -  Mo  -  (M*  +  M  )  o  -  M.6  cos  pt  -  M  ep8  -  0 
x  a  a  q  o  pp 

••  ..  .  (i) 

18  +  palx  +  N08  -  (Nr  -  N|)  8  -  Mfi«  sin  pt  -  N^pa  -  0 
In  the  above  equation,  a  and  8  are  the  angular  deviation  of  the  Model's  axis  of 
symmetry  (or  principal  axis  of  Inertia)  fro*  the  velocity  vector:  a  Measures 
deviation  In  the  vertical  plane  and  8  deviation  in  the  horizontal  plane.  Figure 
6  indicates  the  measurenent  of  these  angles.  Since  the  aaxlaun  excursions  in 
a  and  8  are  limited  to  less  than  7  dagreea,  any  anbigultlea  due  to  connutivity  of 


a  and  8  are  limited  to  less  than  7  dagreea,  any  anbigultlea  due  to  connutl 
angular  rotations  may  be  Ignored. 

Table  1  below  defines  the  various  symbols  appearing  In  equations  (1). 


TABLE  1  ■QQATXGS  (1)  STMOL  DBFXNITIOM 
0 


Cma-Mj/Cqed) 

Cn^/Cqsd) 

CP«-llq/(i^/2V) 

Ch^r/Ciad2/!?) 


Pitching 


eV*pe/<«*4*/lf) 

Gy*4/<5»4> 


sent  derivative 
it  derivative 


Fitch  damping  derivative 
Ten  leg  derivative 


derivative 


Fitch  Manaus  derivative 


Tela 


it  derivative 


HulvaleMt  trla  off-set  angle 


it  of  inertia 


it  of  Inertia 


See  Footnote  3  on  page  19. 


m i 


»  .  *  .  *  »  * «  A  - 


ANGLE  OP 
ATTACK 


FIGURE  6  ORTHOGONAL  CAMERA  COVERAGE  OF  MODEL  MOTION. 
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If  equation  (1)  is  usued  linear  in  angle  of  attack,  a,  end  sideslip,  8  and 
the  spin  rate,  p,  is  constant,  then  various  formalisms  are  available  to  effect  a 
solution.  One  approach  might  be  to  reduce  the  differential  equation  to  an 
algebraic  equation  through  the  Laplace  Tr ana form.  If  the  existence  of  rotational 
symaetry  of  the  model  is  accepted  the  consequence  is  a  set  of  relationships  among 
the  derivatives  of  Table  1  (and  equations  (1»  as: 


Vr 


v  -  V 

The  procedure  taken  here  to  solve  equations  (1)  is  not  to  use  the  Laplace 
Transform  but  rather  to  remain  in  the  time  domain,  first,  equations  (1)  may 
be  combined  by  introducing  the  complex  angle  of  attack,  £,  as 

C  -  ia  +  0  < 

Equations  (1)  can  now  be  expressed  as  a  single  differential  equation  in  the  com¬ 
plex  variable,  £,  as: 

(  +£tp  ^  -  t  ("**  *  *  1  “*  pt)  S  '  0  ‘ 

After  some  protracted  analysis,  the  solution  of  equation  (4)  takes  the  form. 


(A.+iuiJt  (X2+l(o2)t 

€  -  kjC  +  k2e 


+  k3.l»* 


where  the  three  additive  terms  on  the  right  might  be  described  as  vectors  (or 


"arms")  in  the  complex  plane.  The  arms  rotate  at  frequancies  (nutation),  o>2 
(precession)  and  p  (spin  rate).  Because  of  the  real  terms  in  the  exponents  of 
first  two  arms  (nutation  and  precession  respectively)  these  arms  also  diminish  or 
grow  depending  upon  the  sign  of  A^  and  A2<  Since  k^,  k2  and  k^  are  complex,  each 
will  in  general  contain  two  real  constants.  Thus,  writing  the  equivalent  of 
aquation  (5)  in  real  variables,  a  and  8,  we  would  expect  six  unspecified 
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constants.  In  addition,  of  course,  there  are  the  two  damping  terms  X^  and  X^  and 
the  two  frequency  terms  and  u^.  (The  spin  rate,  p,  will  be  treated  as  a 
measurable  quantity.)  Writing  equation  (5)  as  two  real  equations  gives, 

X^t  X  t 

a  =  (A^  cos  w^t  +  sin  ui^t)e  +  cos  o^t  +  B^  sin  u^Oe 

+  (A3  cos  pt  +  B^  sin  pt)  (fca) 

V  V 

8  =  (B^  cos  w^t  -  A^  sin  w^tje  +  (B^  cos  -  A2  s*n  (jJ2t^6 

+  (B^  cos  pt  -  A^  sin  pt)  (6b) 


It  then  appears  that  a  data  reduction  scheme  will  require  the  record  of 
angular  displacements  a  and  8  to  be  "fitted"  in  some  way  to  provide  numerical 
values  of  the  ten  constants:  A^ ,  B^,  k^ ,  B^ ,  A^,  B^,  X^,  ,  w  . 

In  spite  of  the  requirement  of  evaluating  ten  constants  only  four,  X^,  X^ , 
w  ,  w2  have  any  practical  value  in  the  evaluation  of  the  aerodynamic  derivatives. 
The  first  four  constants  A  ,  B  ,  A  ,  B  contain  in  addition  to  the  stability 

A  JL  im 

•  • 

derivatives,  initial  conditions  in  the  variables  a,  a,  8,  3  at  whatever  point  (in 
time)  the  a(t),  8(t)  record  is  begun.  Since  it  is  not  practical  to  evaluate 
precisely  the  above  mentioned  initial  conditions,  the  constants  A^,  B^,  A B2  have 
no  utility  in  establishing  the  stability  derivatives.  For  example,  A^,  may  be 
written  as, 


Ai  ■ 


Eo  '  Vo  -  “2e0  *  *2*3  ~  (P  '  “2)B2](>1  ~  V  (7) 

(wL  -  m2)2  +  (X1  -  X2)2 
The  dependence  upon  initial  conditions  d^,  aQ)  is  clear.  Similar  relationship  might 
be  written  for  B^  A2>  B^  A^  and  B^  are  expressions  for  configurational  asymmetries. 

The  frequency  and  damping  terms  ,  w2 ,  X^,  X2  are  the  constants  of  interest 
since  they  contain  all  the  aerodynamic  stability  derivatives  except  the  a  ymmetry 
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if,  Mgfi.  First  define  the  two  parameters ,  s,  and  t: 


The  parameter,  s,  is  called  the  "gyroscopic  stability  factor"  and  is  a  measure  of 
the  relative  presence  of  aerodynamic  stiffness,  M^/I,  and  gyroscopic  stiffness, 
(plx/21).  For  our  purposes  the  parameter,  t,  will  be  used  to  convey  essentially 
the  same  information,  i.e.  the  relative  Importance  of  spin  or  gyroscopic  effects 
to  aerodynamic  effects.  It  should  be  noted  that  for  statically  stable  missiles 
(center  of  pressure  aft  of  center  of  gravity)  the  stability  derivative,  Mq,  in 
equation  (9)  is  negative.  Thus  t  must  always  be  real.  Further  note  the  asymp¬ 
totic  values  of  t  with  spin  rate,  p. 

11m  r  «*•  0 

p  *►  0 

11m  t  *►  1 


p  «* 

Using  the  parameter,  T,  it  is  easy  to 


show 

(10a) 

(10b) 

(10c) 
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X  — 3.  (1  -  t) 

2  21  ^  ; 


■0) 


Also  It  Is  Important  to  realize  that  the  sign  of  t  Is  that  of  spin  rate,  p. 
Assuming  that  p  Is  positive  It  Is  clear  from  equation  (10a, b)  that  is  positive 
and  b>2  is  negative;  further,  the  magnitude  of  co^  is  greater  than  that  of  u^.  Unde 
conditions  of  positive,  p,  the  greater  frequency  is  and  identified  as  nutation 
and  the  lesser  frequency  is  u>2  and  identified  as  precession.  In  the  data  to  be 
given,  the  spin  rate  is  negative  with  the  result  that  the  nutatlonal  frequency  is 
negative  and  the  precessional  frequency  is  positive.  Regardless  of  the  sign,  the 
frequency  having  the  greater  magnitude  is  the  nutation,  while  the  frequency  having 
the  lesser  magnitude  is  the  precessional. 

Equations  (10)  present  the  four  constants  u^,  oi^,  X. ,  X^  in  terms  of  the 
aerodynamic  derivatives  M  ,  M  and  M  and  the  spin  rate,  p  (through  the  parameter 

°  q  pp 

t).  It  is  relatively  straightforward  to  evaluate  the  terms  M  ,  M  ,  M  .  and  p  (and 

a  q  pp 

t)  in  terms  of  the  four  constants  u^,  «2,  X^,  X^ .  Multiplying  equations  (10a, b) 
together  and  using  equation  (9)  gives: 

M  ■  -Ia),a)_  (11a 

a  12 

Next  adding  equations  (10c, d)  provides: 

Mq  -  I(XX  +  X2)  (lib 

Equations  (I0a,b)  may  then  be  added  to  give: 


p  ■  -J-  (“l  +  w2> 


Next  subtract  equation  (10b)  from  (10a)  and  using  (11c)  above  gives: 

“l  '  “2 

If  equation  (lOd)  is  subtracted  from  (10c)  and  equation  (lib)  is  used,  the 
result  is: 


wwmrriAr rr-w.A 
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M 

X  -  X.  ■  (X  +  X  )t  +  -E2-  (2t) 

1  Z  12  I 

Next  replacing  t  from  equation ■ (lid)  gives  after  aone  manipulation: 


Thus  equations  (lla,b,e)  provide  a  means  for  calculating  the  dimensional  aero¬ 
dynamic  derivatives  Mq(  M^,  from  the  four  constants  u^.  X^,  X^.  In 
addition  it  will  be  noted  that  the  spin  rate,  p,  may  be  calculated  from  the  reduced 
data.  The  conclusion  that  spin  rate  need  not  be  measured  is  a  bit  misleading:  a 
good  estimate  of  spin  rate  is  quite  helpful  In  initiating  parameter  estimation 
during  data  reduction.  After  the  data  reduction  process  has  converged  a  measured 
spin  rate  and  a  computed  spin  rate  from  equation  (11c)  will  Indicate  by  their 
agreement  a  check  on  the  validity  of  the  reduced  data. 

It  will  be  noted  In  equation  (11c)  that  if  the  spin  rate,  p,  is  zero  that 
•  -Wj  •  » 

and  if  p  -  0,  then  t  *  0;  it  follows  from  equations  (10c, d)  that: 


X 


1 


X 


Consequently  equations  (lla,b,e)  become: 
M  -  -Iw2 

ft  • 

M  -  2IX 
9 

M  -  INDETERMINATE 
PP 


(12a) 

(12b) 

(12c) 


Equation  (12c)  indicates  that  in  the  absence  of  spin  rate,  p,  equation  (lie)  takes 
on  an  Indeterminate  form.  It  should  be  appreciated  that  M^,  the  Magnus  derivative, 
is  dependent  upon  configuration  and  not  upon  spin  rate.  However,  in  the  absence 
of  spin  it  is  impossible  to  measure  the  Magnus  derivative. 
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If  the  motion  of  the  model  Is  Anemically  stable,  then  X^,  X,  <  0,  l.e.  both 

the  nutation al  and  p recessional  damping  exponents  must  be  negative.  Hence,  at 

« 

large  values  of  time,  equations  (6)  become. 


aa  -  cos  pt  +  sin  pt 


(13a) 


cos  pt  -  A3  sin  pt 


(13b) 


Assuming  that  the  angles  a m  and  0^  are  small,  it  is  justifiable  to  square  both 
sides  of  the  equations  (13),  add,  and  finally  take  the  root  to  give: 


42+02 


423  +  B2 


where  <*T  is  the  trim  angle  of  attack.  Thus  the  constants  A^  and  will  provide 
an  indication  of  configurational  asymmetries  causing  a  non-zero  trim  angle. 

3.2  DATA  REDUCTION  EQUATIONS.  The  basis  of  the  data  reduction  method  is  a 
variation  of  the  familiar  least  squares  technique.^  In  the  sketch  below  on  the 


te*lT‘1'’ 2ni  “ttlOT  w“ TorkI 


XZ&K"' 
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left  is  shown  a  series  of  angle  of  attack  measurements  and  the  corresponding  value 
of  time.  The  function  a(t)  is  the  best-fitted  function  (of  the  assuned  form)  to  the 
succession  of  data  points  (a^,  t^) .  The  fit  is  "best"  in  the  least-squares  sense. 

In  other  words  for  each  data  point  (a^,  t^)  there  is  a  "residual"  which  for  the 
angle  of  attack  is  designated  by  r^.  Similar  remarks  apply  to  the  angle  of  side¬ 
slip,  8,  shown  to  the  right  in  the  above  sketch.  The  residual  or  the  difference  be¬ 
tween  the  assumed  best-fit  function  and  data  point  at  the  same  time  is  indicated  by 
the  variable,  q^.  Thus  we  may  express  analytically  the  residuals  and  q^  as, 

r^  ■  a(t^;  A^,  B^,  B ^ ,  A^,  B^,  Wj,  X^,  X2)  ”  (15a) 

q±  “  B(tjJ  A^,  B^,  A2,  Bj,  Aj,  Bj,  Wj*  X^,  X2)  —  ^  (15b) 


The  residual  is  shown  as  the  difference  between  the  fitted  curve  and  the  data  point 
at  the  time  t  ■  t^.  The  residuals  and  q^  are  shown  to  depend  upon  the  ten  para¬ 
meters  A^,  .  .  .  Xj*  7116  heart  of  the  least  squares  technique  then  is  the  method 
for  establishing  ten  linearly  independent  equations  for  the  determination  of  these 
ten  constants.  These  equations  are  established  by  minimi  zing  the  sum  of  the  squares 
of  the  residuals  with  respect  to  each  of  the  ten  constants.  Since  there  are  two 
variables  i.e.  a  and  6,  the  quantity  to  be  minimized  is  the  sum  of  each  of  the 
residuals  squared  as. 


(16) 


Then  the  derivative  of  the  above  quantity  is  taken  in  turn  with  respect  to  each  of 
the  ten  constants  A^,  .  .  .  X2;  these  derivatives  are  each  set  equal  to  zero  to  pro¬ 
vide  ten  independent  algebraic  equations  that  can  be  solved. simultaneously  for  each 
of  the  ten  constants.  For  example  for  the  A^  and  B^  constants,  the  derivative 
operation  will  be: 


(17a) 
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t  ?(■!•<)- 


(17b) 


If  the  analytic  expression  for  model  notion  were  linear  in  all  ten  of  the 
constants,  the  ten  equations  which  follow  from  the  differentiation  (see  equations 
(17))  might  be  solved  as  j^iaeav  algebraic  equations.  A  glance  at.  equations  (6) 
shows  that  six  constants  ,  B^,  and  do  vary  linearly.  However,  the 

four  remaining  constants  w^,  w2,  X^,  X2  do  not  vary  linearly:  u^,  u>2  are  arguments 
of  transcendental  functions  and  X^,  X2  are  arguments  of  exponential  functions. 

Thus,  the  direct  method  of  least  squares  must  be  extended  to  an  iterative  method 
known  as  "differential  corrections." 

The  residuals  have  been  designated  in  equation  (15);  the  functions  a(t)  and 
B(t)  given  in  these  residual  expressions  can  be  expanded  in  a  Taylor  series  in  the 
constants  u>2,  X^,  X^ •  The  remaining  six  constants,  A^,  B^,  Aj ,  B2,  A^,  B^ 

could  also  be  evaluated-by  this  iterative  differential  corrections  procedure,  but 
since  the  residuals  vary  linearly  with  these  six  constants,  the  more  direct  least 
squares  procedure  will  be  used  in  their  evaluation.  The  least  squares  method  will 
be  covered  in  more  detail  after  the  differential  corrections  method  has  been 
described.  Replacing  o(t)  and  B(t)  in  equations  (15)  by  a  Taylor  expansion  in  each 
of  the  variables  u>^,  u>2,  X^,  X2  gives: 

r^  *  ®(t^ ;  ,  B^ , . . . ,  ,  w2 ,-  Xx ,  Xj) 

.  3a  .  3a  .  3®  \  j.  I®,  i  .  v 

+  “l  +  KJ  “2  +  3^  X1  *  ITj  X2  *1  (1M 

(J^  *  i  B^i  •  i  •  |  (D^f  ^ i  j  X2) 

+  %  “i + !%■  “j  +  Hj;  li + Ifj  x2  •  8i  (18b> 

where  a  and  3  have  bean  expanded  about  estimated  values  «2,  X^,  The 


(18b) 
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quantities  or  variables  a>^,  <i»2,  X^,  X2  represent  the  Incremental  changes  in  the 
estimated  values  (see  equation  (26)). 

Hew  residuals  and  nay  be  defined  as: 


Ri-o 


(t^»  Aj,n.,  X2)  — 


O-O. 


Ql  -  e(t±S  Aj . X2)  -  0±  “  J  -  &± 


(19e) 

(19b) 


The  residuals  R^  and  are  recognized  is  the  difference  between  die  estimated 
functions  o(t1;  A^,. . . ,  X^)  and  fKt^;  A^,...,  X^)  and  the  date  points  and  8^. 
Equations  (18)  and  (19)  may  be  combined  to  give: 

(20a) 


So  .  So  .  3o  ,  .  lSL-  1  4.  o 
3^  “l  +  3o>2  “2  +  3X1  X1  +  3X2  X2  +  R1 


t 


IS-  w  + _ , 

3u^  3«2  ^ 


30  w  4.  38  -y  +  38__  ^  4.  q 

T  Ai  T  A*»  T  *»4 


ax,  '  1 


ax2  2 


(20b) 


The  residuals  r.  and  q.  may  be  squared  and  summed  to  give  the  quantity 


1 

as  In  identity  (16).  Now  the  derivative  is  taken  of  the  above  sum-of-the-square 
of  the  residuals  with  respect  to  the  four  variables  w^,  &>2’  *2  811(1  tlie8e 

derivatives  set.  In  turn,  equal  to  zero.  For  example  the  equation  which  results 
from  taking  the  derivative  with  respect  to  is: 


30 
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It  should  bs  apprsdstsd  that  the  various  partial  darlvatiuas  in  aquation  (21) 


such  as  srs  evaluated  at  the  estimated  values  «2,  V^,  72  and  hence  are 

constants  for  any  differentiation  with  respect  to  «2,  X^,  X^. 

The  development  of  the  regaining  three  equations  is  obvious.  The  result  is 
that  four  aquations  linear  in  the  variables  «2,  X2  are  available  for  the 
calculation  of  these  variables.  Matrix  algebra  My  be  used  to  advantage  in 
progressing  these  equations  for  solution  on  a  digital  computer.  Let  be  the  Jth 
subscripted  variable  representing  u^,  «2,  X^,  X2  for  the  four  consecutive  sub¬ 
scripts  as 


The  algebraic  equation  given  in  aquation  (22)  and  the  equivalent  equations  for  the 
reaalning  variables  nay  be  represented  by  the  following  natrix  equation: 

c]  ’  “k  <“> 


represents  the  coefficients  of  equation  (23)  and  the  equivalent  equations; 
Rl^  contains  the  sun  of  the  residuals.  For  example  the  4x4  natrix  has  the 
following  elemental  values: 


32 


■;  •'  v*  .  •  »».ri 

_:.-v-V.v:j 
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F11  ■  ^ 

”/  \2  /  \2l 

(— Y  + 1— V 

\  wl/  v  wl/ 

F12’f 

3a  3a  +  38  38  ~j 
3(0^  8(1)2  3w^  3a>2  I 

P13-f 

3a  3a  38  38  ~| 

3«i  3XX  3(o^  3XX 

F14  '  ^ 

'3a  3a  38  38  1 

3W;l  3X2  +  30)2  3X2 

F21  “f 

3a  3a  38  38  1 

3(0  3(o,  3(0.  3(o, 

L  2  1  2  1J 

F22  "  ^ 

(aa  V  ,  4e  VI 

l*v  \»*y  J 

F23  "  ^ 

3a  3a  .  38  38  1 

1“2  ”l  “ll 

P24 

~3a  3a  +  3g_  38  ~| 
3«2  8X2  3»2  3Xj 

F3l’f 

‘3a  3a  38  38  "I 

3X1  3“l  3X1  3“l 

F32*f 

3a  30  .36  38  1 

3XX  3(o2  3Xx  3m2  I 

F33 

F34  ‘  £ 

‘3a  3a  36  38  ~| 

3XX  3X2  +  3XX  3X2  I 

F41-f 

3a  3a  38  38  1 
3X2  3«2  3X2  3Wl  J 

?*2  '  £ 

3a  3a  ,  36  38  "1 

3^  3w2  3X2  3(o2  J 

33 
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_  Vila  3a  .36  36  ~ 

43  ^ 

r«  -?tW + (%H 

The  residual  tan  uy  ba  written  aa 

“i  v + “t] 

“’-?[% •**%*] 

U3  -S[l^  *1  ♦  «l] 

«k-x[%1i+^,i] 


(24a) 


(24b) 


Equations  (26)  say  ba  traatad  aa  a  Matrix  aquation  and  a  solution  for  C.  may  ba 

•J 

Indicated  the  inverse  formalism: 


1  5  c,  -  <V_1  “v '  nik  *H 


Jk  k 


For  most  applications  the  matrix  is  not  singular  so  the  Inverse  exists. 
However  this  matrix  is  not  orthogonal  so  the  inverse  has  been  Indicated  as  Fl^ 
Once  the  variables  u>^,  u>2>  X^,  have  been  obtained  a  now  estimate  on  these 
variables  can  be  made  as 

u>2  +  oij  -*■ 

“2  +  W2  “2 

X^  +  X^  *  X^ 

X2  +  X2  X2 


WPM 
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In  other  words  the  new  estimates  of  the  variables  shown  on  the  right  hand  side  are 
formed  by  adding  algebraically  the  incremental  variables  to  the  present  estimates. 
After  each  iteration  new  estimates  become  available;  these  new  estimates  are  used 
to  reevaluate  the  various  derivatives  appearing  in  the  elements  of  the  F^  and  Rl^ 
matrices  (see  equations  (24)).  After  a  sufficient  number  of  iterations  a  "best 
estimate"  (in  the  least-square  sense)  is  available  for  the  variables  ui^,  u^.  A.,  X^. 
Obviously  selecting  the  number  of  iterations  is  arbitrary.  However,  the  probable 


error  of  fit,  E,  can  be  used  as  an  index  of  convergence.  E  may  be  given  as: 


0.6745 


2  2  ~ll 

r  (g±  + 

~  (2N  -  n) 


where  N  is  the  total  number  of  data  points  to  be  fitted  and  n  is  the  number  of 
parameters.  In  the  case  under  discussion  n  is  obviously  equal  to  four. 

The  derivatives  of  u  and  8  with  respect  to  the  four  parameters  u^,  <i>2»  A^,  A^ 
have  not  been  defined.  However  these  derivatives  are  easy  to  form.  From  equation 
(6)  the  eight  required  derivatives  are: 

■gjj-  *  (t)  (B^  cos  w^t  -  Aj  sin  u^t)e  (28a) 


At 

(t)(B2  cos  «2t  -  A2  sin  «2t)s 


(28b) 


»  Aft  W 

•jj-  •  (t)(A^  cos  w^t  +  B^  sin  u^t)s  1 
da 

-  (t)  (Aj  cos  «i»2t  +  B2  sin  w2t)e 

}«  b 

jJJ-  •  (“t)(B1  sin  «jt  +  A^  cos  Ujt)s 


(28c) 

(28d) 

(28e) 


:  'tti «Ua ft  ■»'*«  *» 


NSWC  TR  81-491 


..  A_l 

■gjj-  "  (-t)(B2  sin  «2t  +  Ajcos  w2t)e  2 

38  V 

■  (t>(B2  cos  u»jt  -  sin  w^t)e 


||-  -  (t)(B2  cos  w2t  -  A2  sin  u>2t)e  2 


(48f) 


(28g) 


(28h) 


Once  a  set  of  values  for  u>2 »  X^,  X2  emerge  from  the  data  reduction 
process ,  It  la  possible  to  calculate  the  aerodynamic  derivatives  and  spin  rate  from 
equations  (lla,b,c,e). 

The  remaining  question  concerns  the  evaluation  of  the  six  constants  A^,  B , , 

Aj,  B2,  Aj,  B3>  which  up  to  now  have  been  assumed  available.  As  pointed  out 
earlier  these  six  constants  must  be  evaluated  to  complete  the  fit  of  equations  (6) 
to  the  data  even  though  such  constants  (with  the  exception  of  and  B^)  contain 
no  new  aerodynamic  data.  It  will  be  noted  In  equation  (6)  that  the  angle  of 
M«ck.  «,  ud  .Id.  .Up,  e,  .r.  Unur  ftmctloo.  of  th.  con.tm.t.  Aj.  »r  *2>  »2> 
Ay  Bjj  hence  the  direct  least  squares  technique  can  be  applied  rather  than  the 
Iterative  differential  correction  method  just  developed  for  the  computation  of 

»  ®2 »  ^2 • 

As  before  (in  equation  (15))  we  may  write  for  the  residuals 

r±  -  o(t)  -  <*t 

q±  -  3(t)  -  81 

where  a(t)  and  $(t)  are  the  functions  given  in  equation  (6)  and  and  0^  are  the 
data  points  at  time  t-t^.  The  functions  a(t)  and  0(t)  are  also  evaluated  at  time 
t»t^.  Of  course  at  this  point  a(t).and  0(t)  cannot  be  numerically  evaluated 
because  A^,  B^,  A2,  iy  A^,  B,  have  not  been  assigned  numerical  values  although  the 
remaining  four  constants,  w, ,  w4,  X, ,  A  have  been  given  estimated  values.  Now  if 
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the  residuals  are  squared  and  added  the  derivative  of  the  resulting  expression  may 

be  taken  in  turn  with  respect  to  each  of  the  unknowns  A^,  B^,  A^,  B^,  ky  B^  to  pro¬ 
vide  six  independent  equations  for  the  evaluation  of  the  unknowns. 

It  will  be  noted  in  equations  (6)  that  the  coefficients  of  the  unknowns  A^»  B^, 
Aj,  B2,  A3,  B3  are  the  derivatives  of  a  and  8  with  respect  to  these  unknowns.  Thus 
equation  (6)  might  be  written  as 

a2  + 


A2  + 


where  the  derivatives  are  given  below  as  : 


3a 

3aT 

v 

cos (u^t)e 

38 

3A1 

3a 

3Bi 

X  t 

sin(a>1t)e 

36 

3BX 

3a 

3A2  " 

At 

cos(u2t)e 

36 

3A2 

3a 

3B2 

At 

sln(a)2t)e 

38 

3B2 

3a 

S' 

cos (pt) 

38 

3A3 

3a 

3B3  “ 

sin(pt/ 

38 

3B3 

cos(h>2t)e 


The  residuals  may  be  summed  over  the  data  points  to  give: 


^  +  "i  ‘  ^  -  °i 


2  +  (B(t)  -  B,]2 


Six  linearly  Independent  algebraic  equations  may  now  be  obtained  by  taking  the 
derivative  of  the  above  residual  sum  with  respect  to  each  of  the  parameters  A. 


i! 
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The  matrix  G  has  36  elements.  Rather  than  writing  all  of  these  elements  a  typical 

K.J 

element  will  be  given;  construction  of  the  remaining  elements  should  be  fairly 

obvious.  First  it  will  be  noted  that  the  column  matrix  has  been  formed  by 

assigning  numbers  1  through  6  in  order  to  the  constants  A^,  B^,  A2,  B^,  A^,  B^; 

then  the  element  G  of  the  matrix  G,  .  is  formed  from  the  derivative  of  a  with 
mn  kj 

respect  to  the  m-th  variable  times  the  derivative  of  a  with  respect  to  the  n-th 
variable  added  to  the  product  of  the  derivatives  of  8  with  respect  to  the  same 
variables.  For  example  the  element  where  m=3  and  n=5,  i.e.  G^,.  >  the  variables  are 


A^  find  A^  (see  column  matrix  in  equation  (35)).  Thus, 


’35 


I 

i 


3a 

38 

33 

3A_ 

3A- 

3A_ 

3AC 

3 

5 

3 

5 

(36) 


Clearly  the  matrix  is  symmetric  (as  is  the  matrix  F^  discussed  earlier  in 
equation  (23) ) .  The  matrix  Qlfc  follows  a  form  identical  to  that  for  Rl^  (equations 
(24b)).  Again  if  k=4,  i.e.  -  B2,  then, 


Ql, 


2 

i 


—  a  +  UL  g 

9B2  i  3B2  Pi 


(37) 


The  elements  of  the  matrix  can  be  computed  by  equation  (35) .  Once  is  avail¬ 
able  then  the  constants  A^,  B^,  A^,  B2>  A^,  B^,  are  also  available.  With  the 
availability  of  these  six  constants  it  is  possible  to  use  the  method  of  differential 
corrections  to  obtain  the  remaining  four  constants  u>^,  o>2>  A^,  The  various 

first  order  partial  derivatives  used  in  forming  the  elements  of  G^  and  Ql^  (see 
equations  (30))  'equired  a  guess  or  estimate  of  the  four  constants  u>^,  b>2.  A,  ,  A^. 
The  method  of  differential  corrections  is  now  followed  to  improve  these  estimates. 
The  iterations  continue  until  convergence  is  sufficient  according  to  some  external 
criterion  such  as  that  stated  in  equation  (27). 

The  main  features  of  the  data  reduction  procedure  might  be  reviewed  in  the 


following  statements 


(a)  Estimate  "reasonable"  values  of  the  constants  a>2 »  . 

(b)  Evaluate  the  derivatives  given  in  equations  (30)  at  each  data 

point. 

(c)  Carry  out  the  products  and  sunmations  to  form  the  elements  of  the 

matrices  and  Ql^. 

(d)  Carry  out  the  matrix  inversion  and  multiplication  as  indicated  in 

equation  (35)  to  solve  for  A^,  B2»  A^,  Bj. 

(e)  Evaluate  the  derivatives  given  in  equations  (28)  at  each  data  point. 

(f)  Carry  out  the  products  and  summations  to  form  the  elements  of  the 

matrices  and  Rl^. 

(g)  Carry  out  the  matrix  inversion  and  multiplication  as  indicated  in 
equation  (25)  to  solve  for  u.*.  Aj* 

(h)  Update  (improve)  the  present  estimates  of  u^,  u>2>  A^,  A^  (equations 

(26)). 

(1)  Return  to  step  (e)  and  repeat  the  process  to  step  (h)  either  a 
preset  number  of  times  or  until  some  convergence  criterion  is  satisfied. 

(j)  Return  to  step  (a)  with  new  estimates  of  o>^,  u>2,  A^t  A^  from  step 
(1)  and  repeat  procedure  until  convergence  satisfies  some  external  criterion. 

In  the  application  of  the  method  discussed  in  this  report  step  (1)  is  omitted. 
In  other  words  after  the  first  update  to  the  parameters  u >2»  A^,  X2  has  been 

completed  an  immediate  return  is  made  to  step  (a). 

3.3  COMPUTER  PROGRAM.  A  computer  program  has  been  written  to  carry  out  the 
data  reduction  procedure  developed  in  the  preceding  section.  This  program  has  been 
listed  in  Table  2.  Where  appropriate  equations  of  section  3.2  will  be  related  to 
the  equivalent  line(s)  in  the  program. 

In  lines  100  through  270  some  descriptive  remarks  are  c.  alned;  most  of 
these  remarks  are  a  brief  restatement  of  the  last  few  paragraphs  of  section  3.2. 
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TABLE  2 

DATA  REDUCTION  PROGRAM 


00100  REN  THIS  PtOGRAM  IS  *  DATA  REBUCTION  PROCEDURE  BASES  UPSN 
OHIO  REN  LEAST  S0UARE8/BIFFERENTIAL  CORRECTIONS  FOR  EVALUATING  THE 
00120  REN  TEN  CONSTANTS  BESCRIBIN6  ROLL-PITCH- YAU  NOTION  HEASUGEB  IN 
00130  REN  THE  UINB  TUNNEL  USINO  THE  THREE-B  SA8  BEARING.  THESE  CONSTANTS 
00140  REN  ARE  (A)  CONSTANTS  OF  INTEGRATION  A1, 01 ,42,12;  (•>  ASYWiCTRY 
00190  REN  CONSTANTS  A3,B3;  CO  NUTAT I 0N/PRCCE88I0N  FREQUENCY 

001 AO  REN  CONSTANTS  U3,U4;  (0)  NUTATION/PRECESSION  BAHPIN6  CONSTANTS 

00170  REN  L3fL4.  II , AND  12  ARE  THE  LONGITUDINAL  AND  TRANSVERSE  N0NENT8 
00180  REN  OF  INERTIA  AND  P  IS  THE  SPIN  RATE.  LI  AND  L2  ARE  THE 
01190  REN  STARTING  AND  ENDING  BATA  POINT  NUNBERS.  VARIABLES 
01200  REN  U3,U4,L3,L4  ARE  FITTED  BY  DIFFERENTIAL  CORRECTIONS. 

00210  REN  PROCEDURE  IS  TO  ESTINATE  THE  VARIABLES  U3,U4,L3,L4  AND 
00220  REN  CALCULATE  Al ,11 ,42,02,43,03  FRON  STRAIGHT  LEAST  SQUARES  ; 

00230  REN  USING  ESTINATES  ON  01,01,02,02,03,03  FRON  LEAST  80IIARE8, 

00240  REN  DIFFERENTIAL  CORRECTIONS  ARE  USED  TO  INPROVE  ESTIHATE8 
00250  REN  OF  U3,U4,L3,L4.  THIS  PROCEDURE  HAY  BE  REPEATED  AN  ARBITRARY 
00240  REN  NUNBER  OF  ITERATIONS.  A  PRINT  OUT  OF  THE  STANDARD  DEVIATIONS, 

00270  REN  18  USED  A  A  NEASURE  OF  THE  DEGREE  OF  FIT. 

00280  II  «  .02474 
00290  12  -  0.1925 
00300  BASE  1 

00310  BIN  A(1000),B(1000) , A0( 1 000), B0( 1 000), T( 1000) 

00320  BIN  6(4,41,61 (4,4>,F(4,4),Ft (4,4) 

00330  DIN  R1(4), 01(4), 0(4), R(IOOO) ,0(1000) ,C(4) 

00340  FILE  G1  ■  "TAPE1" 

00350  INPUT  B1,R9,N 

00340  PRINT  "ENTER;  L1«  STARTING  LINE  NO  AND  L2  -  ENDING  LINE  NO.* 

00370  INPUT  L1,L2 

00380  8  ■  LI  -1 

00390  IF  S<»0  THEN  00430 

00400  FOR  I  •  1  TO  8 

00410  INPUT  I1,B,D,D 

00420  NEXT  I 

00430  N  «  L2-LH1 

00440  FOR  I  •  1  TO  N 

00450  INPUT  Nl ,T(I),A(I),B(I) 

00440  NEXT  I 

00470  PRINT  "STARTING  AND  ENDING  VALUES" 

00400  PRINT  T(1 ),A(t ),B(1 ) 

00490  PRINT  T(N),A(N),B(N) 

00500  U3  *  10.42 
00510  U4  -  (-1 )*13.4S5 
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TABLE  2  (CONTINUED) 

DATA  REDUCTION  PROGRAM 

0*520  L3  -  <-!)*. 12 
0*530  L4  ■  <-t>*.0475 
0*540  P  -  (I2/I1)*(U3+U4) 

0*550  FOP  K  ■  t  TO  15 

0*540  BEF  FNA(T)  >  C0S<U3*T)*EXP<L3*T) 

0*570  OEF  FNB(T)  -  8IN<M3*T>*EXP<13*T> 

00580  OEF  FNC(T)  •  C08(U4«T)*EXP(L4*T) 

00590  OEF  FND(T)  *  SIN(U4*T)*EXP<14*T) 

00400  OEF  FNE(T)  ■  C08(P*T) 

0*410  OEF  FNF(T)  *  SIN<P*T) 

0*420  OEF  FN6(T)  »  <04*C08<U3*T)-A4*SIN(U3*T>)*T*EXP<13*T) 
0*430  OEF  FNN(T)  >  (85*C08(U4*T)-A5*SIN<IM*T))*T*EXP<L4*T) 
0*440  OEF  FNKT)  >  <A4*C0S(U3*TK84*SIN<M3*T)>*T*EXP<L3*T) 
0*450  OEF  FNJ(T)  *  <A5«C0S(U4»T)»B5*SItt<U4*T>>*T*EXP<L4*T> 
0*440  OEF  FNM(T)  -  (-1 )*SIX(U3*T)*EXP(L3*T ) 

0*470  OEF  FNN(T)  -  C0S<M3*T)*EXP(L3*T) 

0*480  OEF  FNP(T)  >  (-1 )*8IH(W4*T)*EXP(L4*T) 

0*490  OEF  FNO(T)  >  C08(U4*T)*EXP(L4*T) 

0*700  OEF  FNR(T)  -  (-1 >*SIN(P*T> 

0*710  OEF  FNS(T)  -  COS(P*T) 

0*720  OEF  FNT(T)  -  (-T)*(B4*8I8<03*T)+M*C08(U3*T))*EXP(L3«T) 

0*730  OEF  FNU(T)  *  (-T)*<85*SI*<U4*T)+A5*C0S<M4*T) )*EXP(L4«T) 

0*740  OEF  FNV(T)  *  <T)«(84«C0B<tf3«T)-M»8IN(N3*T))*EXP(L3*T) 

0*750  OEF  FNH(T)  -  (T)*(03*C08(U4*T)-A5*8I8(«4<*T))*EXPa4*T) 

0*740  FOP  I  ■  1  TO  4 

0*770  FOR  J  *  1  TO  4 

0*700  0(1, J)  ■  0 

0*790  NEXT  J 

00000  NEXT  I 

00810  FOP  I  -  1  TO  4 

00020  OKI)  -  0 

00030  NEXT  1 

00040  FOP  I  •  I  TO  N 

00050  N  •  T(I) 

00040  0(1,1)  -  FN4(H)*FN4(H)+FNN(H)*FNH(N)+G( 1, 1) 

00070  0(1,2)  -  FN0<N)*FNB<NKFNH<H)*FNN<H)-»6<1,2> 

OOOOO  0(1,3)  -  FN0<N)»FNC(H)*FNN<M)«FNP(H)*B<t,3> 

00090  0(1,4)  •  FN4(N)PFN0(H)*FNN(H)*FN0(H)*O(1,4) 

00900  0(1,5)  •  FN4(H)*FNE(H)*FN«(M)*FNR(H)48(1,5) 

00910  OK ,4)  •  FN4(H)*FNF(H)*FNH(H)*FN8(H)40(1,4) 

00920  0(2,1)  •  FNB(N>*FNA(HKFNN(H)*FNH(HK8(2,1) 

00930  0(2,2)  -  FN8(M)«FN8<H»FNN(H)'FNN(H)+8<2,2) 


42 
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TABLE  2  (CONTINUED) 
DATA  REDUCTION  PROGRAM 


01940  8(2,3)  >  FNB(H)*FNC(H)+FNN(H)*FNP(N)tO(2,3) 
00950  8(2,4)  -  FNB(H)*FNB(H)+FNN(H)*FNQ(N)+G(2,4) 
00900  8(2,5)  -  FNI(N)*FNE(N)+FNN(1I)*FNR(H)*0(2,5) 
00970  6(2,0)  -  FNI(H)*FNF(N)+FNN(N)*FNS(H)+6(2,0) 
00980  8(3,1)  -  FNC(H)*FNA(H)+FNP(H)*FNN(H)*G(3, 1) 
00990  8(3,2)  -  FNC(H)*FNB(H)+FNP(H)*FMN(N)*8(3,2) 
01000  6(3,3)  *  FNC(N)*FNC(r*FNP(N)*FNP(H)+8(3,3) 
01010  8(3,4)  -  FNC(H)*FNB(H)+FNP(N)*FNQ(H)+8(3,4) 
0(020  6(3,5)  -  FNC(N)*FNE(N)+FNP(»:)*FNR(N>+G(3,3) 
01030  8(3,0)  -  FNC(M)tFNF(H)+FNP(N)*FN8(H)+8(3vA) 
01040  6(4,1)  •  FNI(H)*FNA(H)+FN8(N)*FNK(H)+0(4,1> 
01050  8(4,2)  >  FNB(H)*FNB(H)+FN8(H)*FNN(N)+G<4,2> 
01000  8(4,3)  >  FNI(N)*FNC(N)+FN8(N)*FNP(H)+8(4,3) 
01070  6(4,4)  ■  FNB(H)*FNB(H)+FN8(H)*FNQ(H)*8(4,4) 
0(080  6(4,5)  >  FNI(N)*FNE(H)*FN8(N)*FNR(H)*8(4,5) 
0(090  6(4,0)  *  FNB(H)*FNF(H)*FN8(H)*FN8(H)+G(4,4) 
0(100  6(5,1)  *  FEE (H)*FNA(H)+FNN(H)*FNH(H) +8(5,1) 
oiiio  8(5,2)  ■  FWE(ir)»FirB(H)»FinmnynHi(nr+gc572i 
01120  6(5,3)  -  FEE(H)*FEC(H)+FEt(H)*FNP(N)+G(5,3) 
01130  8(5,4)  •  FNE(N)*FNI(H)+FNft(N)*FNQ(H)+6(S,4) 
01140  8(5,5)  -  FNE(H)*FNE(H)«FNR(N)*FNR(H)+8(5,5) 
01150  6(5,0)  -  FNE<H)»FNF(H)+FNR(H)»FN6(N)+G(5,0) 
01100  6(0,1)  ■  FEF(N)*FEA(H)+FE8(N)*FEH(N)+8(0, 1 ) 
01170  8(0,2)  ■  FNF(N)*FNB(H)+FN8(H)*FNN(H)+6(0,2) 
01180  6(0,3)  ■  FEF(E)*FEC(N)+FE8(N)*FNP(N)+G(0,3) 
01190  6(0,4)  -  FNF(H)*FEI(N)+FE8(N)*FEQ(H)+B(0,4) 
01200  6(0,5)  ■  FNF(H)tFNE(H)+FNS(N)*FNR(H)+8(0,5) 
01210  6(0,0)  -  FEF(N)*FEF(H)+FE8(E)*FN8(E)+G(0,6) 
01220  61(1)  -  FNA<H)*A(I)+FNH(N)*B(1)+61(1) 

01230  81(2)  •  FN0(N)*A(I)+FNN(H)*»(I>+G1(2> 

01240  81(3)  -  FNC(N)*A(I)+FNP(H)t>tI)481(3) 

01250  81(4)  -  FNI(N)*A(I)+FN8(N)*0(I)+61(4) 

0(200  81(5)  ■  FNE(N)*A( ! )+FNE(H)*t(I )+81 (5) 

01270  81(0)  -  FNF(N)*A<I)+FNS(N)*0(I)+01(0) 

01280  NEXT  I 
01290  NAT  81  -  189(B) 

01300  NAT  8  •  01*81 
01310  A4  -  0(1) 

01320  14  •  1(2) 

01330  AS  -  1(3) 

01340  05  •  0(4) 

01350  AO  -  0(3) 
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TABLE  2  (CONTINUED) 
DATA  REDUCTION  PROGRAM 


01340 

01370 

01380 

01300 

01400 

01410 

01420 

01430 

01440 

01490 

01440 

01470 

01480 

01400 

01900 

01910 

01520 

01930 

01940 

01990 

01540 

01570 

01580 

01500 

01400 

01410 

01420 

01430 

01440 

01490 

01440 

01470 

01480 

01400 

01700 

01710 

01720 

01730 

01740 

01790 

01740 

01770 


84  -  1(4) 

FOR  I  -  1  TO  N 

T  *  T (I ) 

40(2)  -  ( A4*C0S <U3*T ) +14*8 IN (U3*T) >*EXP< L3*T ) 

40(1)  -  40(I)*(45*C0S(U4*mi5*SIN(U4*r»*EXP(L4*T> 
40(1)  •  40(1)  +44*C08( P*T ) ♦14*8111  ( P*T ) 

10(2)-  ( 04*008 ( U3*T ) -44*8 IN( W3*T ) ) *EXP (L3*T) 

80(2)  -  80( I ) + ( 85*C0S ( IM*T ) -45*SIN( U4*T ) )*EXP (L4*T) 
80(1)  -  80(I)*84*C0S<P*T)-44*81N(P*T> 

NEXT  2 

NOT  X  -  40-4 
NOT  0  *  80-B 
FOR  2  -1  TO  4 
FOR  J  ■  1  TO  4 
F(ItJ)  •  0 
NEXT  J 
NEXT  2 

FOR  I  •  1  TO  4 
R1(2)  -  0 
NEXT  I 
ZO  -  0 

FOR  2  •  1  TO  N 
N  -T(2) 

F(1V1)  •  FN8(N)*FNQ(H)*FNT(H)*FNT(H)*F(1,1) 

F(1f2)  -  FN8(H)*FNH(H)+FNT(H)*FNU(H)-fF(1 ,2) 

F(1r3)  »  FNB(H)*FNI(H)*FNT(H)*FNV(H)*F(1,3) 

F(1 ,4)  -  FN6(H)*FNJ(H)+FNT(H)*FNU(H)+F(1 ,4) 

F(2,1)  -  FNH(H)*FN8(H)+FNU(H)*FNT(H)+F (2,1 ) 

F( 2,2)  *  FNN(H)*FNH(H)+FNU(H)*FNU(H)*F(2t2) 

F(2,3)  •  FNH(H)*FNI(H)+FNU(H)*FNV(H)+F(2,3) 

F(2|4)  -  FNH(H)*FNJ(H)+FNU(N)*FNU(H)+F(2y4) 

F(3»1)  -  FNI(H)*FN6(N)«FNV<H)*FNT(H)*F(3,1) 

F(3f2)  •  FNI(N)*FNN(N)*FNV(H)*FNU(H)*F(3t2) 

F(3,3)  *  FNI(K)»FNI(N)«FNV(N)*FNV(H)+F(3,3) 

F(3,4)  »  FNI(H)*FNJ(H)+FNV(H)*FN«(N)*F(3,4) 

F(4,1)  -  FNJ(N)*FN0(N)*FNU(H)«FNT(M)*F(4,1) 

F(4,2)  •  FNJ(H)*FNN(N)*FNM(N)*FNU(H)+F(4,2) 

F(4,3)  •  FNJ(H)*FNI(N)+FNU(H)tFNV(H)+f (4,3) 

F(4»4)  •  FNJ(H)*FNJ(H)*FN«(H)*FNU(H)*F(4»4) 

R1(1)  -  (-1 )*(FN8(H)*R(I)*FNT(H)*8(I))+R1(1) 

R1(2)  «  (-1)*(FNM(H)*R(I)*FNU(H)*0(2))fR1 (2) 

R1(3)  -  (-1)*(FN2(N)*R(I)*FNV(H)«9(2))+R1(1) 


44 
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TABLE  2  (CONTINUED) 
DATA  REDUCTION  PROGRAM 


*4 

*-  I 


01780  R1<4>  *  (-1 )*(FNJ(H)*R(I)*FNU(H)*Q(I))+R1 (4) 

01790  Z8  •  R(I)*R(I)+B(I)*fl(I)*Z8 

01800  NEXT  I 

01810  NOT  FI  «  INV(F) 

01820  NOT  C  >  F1*R1 
01830  U3  -  U3*C(  1 ) 

01840  114  >  M4*C(2) 

018S0  L3  *  L3+C(3) 

01860  L4  >  L4+C(4) 

01870  E  *  0.4745*8QR(Z8/(2*N-4» 

01880  IF  K>1  THEN  01930 
01890  PRINT  "RUN  NUNIER  *MfR9 
01900  PRINT 

01910  PRINT  "ST AID. DEV. -E  "f "NU-FRE8" , "PREC-FREO" t "NU-IOHP" f "PREC-IANP" 
01920  PRINT 

01930  PRINT  E,U3,U4,L3,L4 
01940  P  «  (12/11 )*(U3+W4) 

01930  NEXT  K 
01960  PRINT 

01970  PRINT  *81V*1  VR2Vi2,,l,,83B 
01980  PRINT  M,84,A5tl5,A6 
01990  PRINT 

02000  PRINT  "83"  f  "111  **,  “U2"  »  *L1  "y  **12" 

02010  PRINT  86^3,84,13,14 
02020  EN8 
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In  lines  280  and  290  the  transverse  and  longitudinal  moments  of  inertia  of  the  model 
are  given.  Line  300  may  not  be  necessary  with  all  computer  systems;  here  the  term 
"BASE  I"  indicates  that  all  subscripted  variables  begin  with  "one"  rather  than 
"zero."  Some  BASIC  compilers  assume,  in  default,  that  a  matrix  dimensioned  D(4) 
has  five  elements.  The  statement  of  line  300  means  that  a  dimension  of  D(4)  has 
four  elements. 

Various  matrices  are  dimensioned  in  lines  310  to  340.  The  matrices  A,  B,  AO, 

BO  and  T  have  been  dimensioned  in  excess  of  any  encountered  data  length.  However, 
all  five  column  matrices  should  have  the  same  dimension  number.  This  remark  also 


applies  to  column  matrices  R  and  Q  in  line  330.  All  other  matrices  are  dimensioned 
to  the  exact  value  needed  in  the  program. 

Lines  340  through  490  are  for  data  input.  All  material  in  these  lines  can  be 
omitted  provided  that  some  other  way  is  chosen  to  input  data.  Data  must  be  input 
as  three  column  matrices:  angle  of  attack,  a(tj)  “  A(I) ,  angle  of  side  slip, 

BCtj)  ■  B(I)  and  time  t^  -  T(I),  where  I  is  data  count.  The  data  triplet  A(I), 

B(I),  T(I)  need  not  be  at  equal  time  intervals,  but  the  triplet  must  be  complete. 

In  the  illustrated  arrangement,  data  is  available  from  a  file.  The  symbol  R9  in 
line  350  is  used  to  designate  the  file,  here  taken  as  a  programmer's  identification 
number.  The  symbol  tt,  also  in  the  same  line,  indicates  the  number  of  data  point, 
i.e.  how" many  triplets  a(t),  B(t),  t  are  in  the  list.  In  line  370,  an  input  state¬ 
ment  permits  restricted  input  from  the  data  file.  Thus  LI  is  the  data  point 

number  at  which  the  file  begins  and  L2  is  the  data  point  number  at  which  the  file 
ends.  To  enter  the  total  data  file,  Ll  ■  1,  L2  -  N  (numerical  equal  to  total  number 
of  data  points). 

In  Unas  500  through  530  initial  estimates  of  the  four  constants  -  W3, 

«2  ■  W4,  X^  ■  L3  and  X2  ■  L4,  according  to  step  (a)  in  the  summary  given  at  the  end 


of  section  3.2.  An  estimated  value  of  spin  rate,  p,  is  given  in  line  540,  which  is 
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based  upon  the  estimates  of  u>^  and  «2  and  equation  (11c) . 

Line  550  sets  the  number  of  Iterations  that  the  program  will  undergo.  Obvi¬ 
ously  this  number  is  arbitrary;  a  conditional  transfer  might  be  entered  near  the 
end  of  the  program  based  upon  some  selected  value  of  E;  however,  this  was  not  done. 

The  functions  defined  in  lines  560  through  750  are  easily  related  to  the 
twelve  derivatives  given  in  equation  (30) .  Lines  860  through  1270  define  the  36 
elements  of  the  G.  ,  matrix  and  the  six  elements  of  the  01 .  matrix  of  equations  (34) . 

kj  j 

Both  of  these  matrices  are  initialized  in  lines  760  through  830.  The  solution  of 
the  matrix  equation  (equation  (35))  is  given  in  line  1300.  The  elements,  i.e. 

A^,  B  ,  A2,  B2,  A3>  B^  are  obtained  in  lines  1310  through  1360;  A^,  B  ,  A2,  B  ., ,  A, , 
B3  are  equal  in  turn  to  A4,  B4,  A5,  B5,  A6,  B6. 

In  lines  1390  through  1440  the  initial  estimates  on  a(t)  and  8(t)  are  given. 
According  to  equations  (18)  these  initial  estimates  are  represented  by  a(t^;  A^, 

...,  X2)  and  8(ti;  A^,...,  X2)  and  are  formed  by  using  the  estimated  parameters  in 
equations  (6).  The  residuals  and  of  equations  (19)  are  set  in  lines  1460  and 
1470.  The  elements  of  the  four  by  four  matrix,  F^ ,  of  equation  (23)  are  calculated 
in  lines  1590  through  lines  1780.  The  various  functions  used  in  calculating  the 
elements  of  are  defined  in  lines  1620  through  1650  and  lines  1720  through  1750. 
The  Rl^  column  matrix,  defined  in  equations  (24b)  is  calculated  in  lines  1750 
through  1780.  The  matrix  inversion  indicated  in  equation  (25)  is  carried  out  in 
line  1810.  The  elements  of  matrix  are  calculated  in  line  1820.  The  Iterations 
of  equation  (26)  is  given  in  lines  1830  through  1860  where  ui^,  o>2,  X^,  X2  are 
represented  in  turn  by  W4,  W5,  L4,  L5.  The  probable  error  of  fit  expressed  in 
equation  (29)  is  calculated  in  line  1930. 

The  iterative  process  is  repeated  by  Indexing  K  in  line  1950.  The  index 
variable,  K,  is  set  at  20  in  line  550.  As  pointed  out  earlier  the  number  of  itera¬ 
tions  are  arbitrary  and  might  in  some  applications  be  keyed  to  the  probable  error 


NSWCTR  81-491 


TABLE  3 

STABILITY  DERIVATIVE  PROGRAM. 


00100  REM  THIS  PRIOR AM  CALCULATES  THE  STATIC  STABILITY, 

00110  REN  DYNAMIC  STABILITY  AN!  NA0NU8  DERIVATIVES, 

00120  REH  SPIN  RATE  AND  CENTER  IF  PRESSURE  FR!H  THE  MUTATIONAL  AMI 
00130  REN  PRECESSIINAL  FREQUENCIES  AND  EXPONENTIAL  DAHPINO. 

00140  REN  UIND  TUNNEL  NACH  NUMBER  (H),  TOTAL  PRESSURE  <PO),PSF, 
OOISO  REN  TOTAL  TEMPERATURE  (TO),  DEO-F,  NUTATI0NAL/PRECE88I0NAL 
00140  REN  FRE0UENCIES/DANP1N6  U1,U2,L1,L2  ARE  ENTERED.  D,I,I1 
00170  REN  ARE  REFERENCE  LENSTH  (DIAMETER), FT,  AND  TRANSVERSE  AND 
00180  REN  AXIAL  MOMENTS  OF  INERTIA.  N  IS  THE  NORNAL 
00190  REH  FORCE  DERIVATIVE  (PER  RADIAN). 

00200  HI  -  14.2148 

00210  U2  >  (>t)«10.3401 

00220  LI  >  (-1 )*2. 15441E-2 

00230  L2  *  (-1 )*2.53824E-2 

00240  N  -  1.35 

00250  B  >  .5045 

00240  I  -  .02474 

00270  II  «  .1925 

00200  N  «  7.95 

00290  PO  »  1.0144E5 

00300  TO  •  955 

00310  S  •  3.14159*(D*D)/4 

00320  TO  »  TO+459.4 

00330  0  »  < ( 1 ♦<N42)/5) ♦( (-1)47/2) )*0.7* (H92)*P0 
00340  V  »  49.0t43*S0R(T0) 

00350  V  «  V4N*((H(Nf2)/5)f((-1)/2>) 

00340  N1  -  IUU14U2 
00370  N2  «  I1*(L1+L2) 

00300  N3  •  (-1)*I*((L14U2+L24U1)/(U1*U2)) 

00390  P  -  (11/1)4(01402) 

00400  Cl  ■  N1/(0484B) 

00410  C2  ■  24N24V/(0484|4B) 

00420  C3  ■  24N34V/(04S4|4») 

00430  X  •  (M)4C14.3059)/N 
00440  PI  ■  P*D/(2*V) 

00450  PRINT  "STATIC", "DYNAMIC", "NAONUS", "SPIN", -C.P." 

00440  PRINT  " CNA " , "CNO" ,"CNPA","P( RAD/SIC ) " , "X( M1Y -LENOTH ) ” 

00470  PRINT 

00480  PRINT  Cl ,C2,C3,P,X 
00490  END 


tA.-.-S 
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of  fit,  E. 


The  computer  program  just  discussed  provides  the  nutational  and  precesslonal 

frequencies  and  and  the  nutational  and  precesslonal  exponents,  X^  and  X^. 

However  these  four  quantities  are  only  intermediate  values;  final  data  must  be  In 

the  form  of  the  aerodynamic  derivatives:  the  pitching  moment  derivative, 

Cm  ,  the  pitch  damping  derivative,  Cm  ,  the  Magnus  derivative,  ;  and  the  spin 
a  q  pa 

rate,  p.  A  second  program  has  been  written  to  convert  the  nutational  and  pre- 

cessional  frequencies  and  damping  exponents,  along  with  relevant  wind  tunnel  flow 
measurements  into  these  derivatives.  This  program  is  listed  in  Table  3. 

In  Table  3  the  frequencies  and  damping  exponents  are  entered  in  lines  190 
through  220  with  «  ,  u>2,  X^,  X^  equal  in  turn  to  Wl,  W2,  LI,  L2.  Model  diameter 
is  entered  in  line  230.  Axial  and  transverse  moments  of  inertia  are  entered 
respectively  in  lines  240  and  250.  Mach  number,  total  pressure  and  total  tempera¬ 
ture  are  entered  in  lines  260,  270  and  280  respectively.  The  derivatives  C  , 

®a 

C  ,  C  ,  p  and  pd/2V  printed  out  in  order  in  line  450. 
n  m 
q  pa 

In  Table  4  to  follow  the  symbols  used  in  this  analysis  are  given,  both  the 
symbol  as  it  appears  in  the  analytic  expressions  and  as  it  appears  in  the  computer 
program.  The  number  of  the  first  equation  in  which  the  symbol  appears  and  the 
first  line  in  which  it  appears  in  the  program  are  also  given.  Line  numbers  indi¬ 
cated  with  an  asterisk  refer  to  the  computer  program  given  in  Table  3;  lines 
without  such  a  designation  refer  to  the  data  reduction  program  given  in  Table  2. 

In  the  next  section  some  sample  wind  tunnel  data  is  presented.  These  data 
will  indicate  the  utility  of  the  analytic  methods  just  developed  in  reducing  data 
from  the  angular  motion  of  gas-bearlng-supported  models. 
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TABLE  4 


NALYTIC 

COMPUTER 

FIRST 

FIRST 

SYMBOL 

SYMBOL 

EQUATION 

LINE 

A1 

A4 

(6a, b) 

1310 

A2 

A5 

(6a, b) 

1330 

A3 

A6 

(6a, b) 

1350 

B1 

B4 

(6a ,b) 

1320 

B2 

B5 

(6a, b) 

1340 

B3 

B6 

(6a, b) 

1360 

CJ 

C 

(23) 

1820 

C 

Cl 

380* 

“a 

C 

C2 

390* 

c 

C3 

400* 

-PB 

d,D 

D 

“ 

230* 

dj 

D(J) 

(35) 

1300 

E 

E 

(27) 

1870 

F(J,K) 

(23) 

1590 

GJk 

G(J,K) 

(34) 

860 

i.iy 

12 

(4) 

280 

*x 

11 

(4) 

290 

M 

M 

— 

260* 

n 

N 

— 

350 

— 

1 

— 

400 

— 

J 

— 

770 

K 

550 

SYMBOL  DEFINITION 

DEFINITION 

Initial  Condition  (sym.)  in  a 
Initial  Condition  (sym.)  in  a 
Initial  Condition  (asym)  in  a 
Initial  Condition  (sym)  in  3 
Initial  Condition  (sym)  in  $ 

Initial  Condition  (asym)  in  8 
Frequency/damping  matrix 
Pitching  Moment  derivative 
Pitch  damping  derivative 
Magnus  derivative 
Body  Diameter 
Initial  Condition  Matrix 
Probable  Error  of  fit 
Normal  matrix-differential  corrections 
Normal  matrix-least  squares 
Transverse  moment  of  Inertia 
Axial  moment  of  Inertia 
Mach  number 
Number  of  data  points 
Index  variable 
Index  variable 
Index  variable 


540 


Spin  rate 


P 

pd/2V 


PI 


(11c) 


410*  Reduced  spin  rate 


TABLE  4  (Cont.) 


ANALYTIC 

COMPUTER 

FIRST 

FIRST 

SYMBOL 

SYMBOL 

EQUATION 

LINE 

P 

PO 

270* 

o 

q 

Q 

— 

310* 

qi 

— 

(15b) 

— 

Qi 

Q(l) 

(19b) 

1470 

Q1k 

Q1(K) 

(34) 

1220 

ri 

— 

(15a) 

— 

Ri 

R(I) 

(19a) 

1460 

“k 

R1(K) 

(23) 

1750 

s 

S 

— 

290* 

t 

T(I) 

(4) 

310 

T 

TO 

300* 

o 

V 

V 

— 

320* 

a 

A 

(6a) 

310 

a 

AO 

(19a) 

1390 

a 

00 

“T 

— 

— 

— 

0 

B 

(6b) 

310 

B 

BO 

(19b) 

1420 

8 

MW 

MM 

00 

X 

— 

(12b) 

— 

h 

L4 

(10c) 

520 

h 

C(3) 

(26) 

1850 

X2 

L5 

(10d) 

530 

I2 

C(4) 

(26) 

1860 

DEFINITION 

Total  pressure 

Dynamic  pressure 

Residual  -  least  squares 

Residual  -  differential  corrections 

Residual  matrix  -  least  squares 

Residual  -  least  squares 

Residual  -  differential  corrections 

Residual  matrix  -  diff.  corrections 

Body  cross-sectional  area 

Time 

Total  temperature 

Airspeed 

Angle  of  attack 

Estimate  of  angle  of  attack 

Steady  state  angle  of  attack 

Magnitude  of  trim  angle 

Angle  of  side  slip 

Estimated  angle  of  side  slip 

Steady  state  angle  of  side  slip 

Zero-spin  damping  exponent 

Nutatlonal  damping  exponent 

Estimate  of  Nutatlonal  damp,  exponent 

Processional  damping  exponent 

Estimate  of  Processional  damp,  exponent 
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TABLE  4 

(Cont . ) 

ANALYTIC 

SYMBOL 

COMPUTER 

SYMBOL 

FIRST 

EQUATION 

FIRST 

LINE 

DEFINITION 

5 

— 

(3) 

— 

Complex  angle  of  attack 

T 

— 

(9) 

— 

Inertia/Aerodynamic  Coupling  parameter 

(0 

— 

(12a) 

— 

Zero-spin  oscillation  frequency 

"l 

W3 

(10a) 

500 

Nutational  frequency 

“l 

C(l) 

(26) 

1830 

Estimate  of  nutational  frequency 

“2 

W4 

(10b) 

510 

Precessional  frequency 

“2 

C(2) 

(26) 

1840 

Estimate  of  nutational  frequency 

K 
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FIGURE  7  EXTERNAL  DIMENSIONS  OF  WIND  TUNNEL  MODEL. 
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4.$  WIMP  TUNNEL  TEST 

A  series  of  wind  tunnel  tests  were  carried  out  using  the  sphere  cone  model 
depicted  in  Figure  1.  The  testing  was  all  done  in  the  Naval  Surface  Weapon 
Center's  Hypersonic  Wind  Tunnel.  In  these  tests  the  model  had  both  a  non-ablating 
and  ablating  spherical  nose.  A  line  drawing  of  the  model  is  given  in  Figure  7. 

It  will  be  noted  that  the  hemispherical  nose  has  a  bluntness  of  0.223  (nose  radius 
divided  by  cone  base  radius).  The  model  oscillates  about  an  axis  that  is  0.595 
body  lengths  aft  of  the  unablated  nose. 

Wind  tunnel  tunnel  testing  of  the  sphere-cone  model  is  a  fairly  straight¬ 
forward  procedure.  The  model  is  mounted  to  the  gas  bearing  assembly  shown  in 
Figure  2 .  The  configuration  as  tested  is  shown  fixed  to  the  axial  sting  support 
in  figure  1 . 

With  the  model  mounted  in  the  wind  tunnel,  the  spin  turbine  is  activated  by 
admitting  air  into  the  turbine  conduit  (see  Figure  5b) .  A  tachometer  reading  is 
used  to  monitor  the  spin  rate.  Minor  adjustments  are  made  to  turbine  air  pressure 
to  maintain  spin  rate  after  wind  tunnel  airflow  is  established.  The  angular 
transducers  are  continuously  recorded  and  the  orthogonal  cameras  are  put  into 
operation  just  prior  to  initiation  of  the  kicking  jet. 

Because  it  is  not  possible  to  visually  Inspect  the  model  during  the  period  of 
wind  tunnel  operation,  it  was  found  convenient  to  enter  the  outputs  from  the 
horizontal-6  and  vertical-a  angular  transducers  into  the  XY  jacks  of  an  oscillo¬ 
scope.  The  display  is  a  rough  portrait  of  the  a-6  motion  of  the  model.  This 
practice  enables  the  project  engineer  to  observe  any  angular  instabilities  as  well 
as  the  more  likely  damped  motion.  When  the  model  has  damped  to  a  sufficiently 
SkJAll  angle  or  to  an  obvious  trim  condition,  the  kicking  Jet  is  agsin  used  to 
activate  angular  motion. 
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FIGURE  8  REPRESENTATIVE  MODEL  PHOTOGRAPHIC  VIEWS  APPROXIMATELY  1  SECOND  APART. 
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It  «u  acatloned  earlier  that  while  the  fiber  optic  eysten  functioned,  the 
performance  was  not  adequate  for  data  reduction.  Consequently  the  photographic 
record  froa  the  orthogonal  camera  Installation  was  used.  Figure  6  Indicates  how 
cameras  were  positioned.  One  camera  obtains  angular  motion  of  the  model  in  the 
vertical  plane,  the  other,  model  motion  in  the  horizontal  plane.  After  film 
reading  and  time-frame-count  correlation  the  two  records  a^(t)  and  8^(t)  are 
available  for  the  data  reduction  procedure  described  In  the  previous  section. 

4.1  NOH-ABLATING  HOSE  (MAN)  WIMP  TUNNEL  TESTS.  In  the  non-ablating  nose 
(NAN)  tests  the  7-degree  conical  model  was  tested  with  a  steel  hemispherical  nose 
In  place.  Hie  external  geometry  of  the  wind  tunnel  model  Is  shown  In  Figure  7. 

Hie  testing  procedure  was  Identical  to  that  described  in  section  4.0.  Primary 
data  In  this  test  was  over  900  frames  of  16mm  film  of  model  motion  measured  at 
each  of  the  two  orthogonal  cameras,  typical  frames  of  the  horizon al  (g)  and 
vertical  (a)  cameras  are  presented  In  Figure  8  .  Hie  total  of  more  than  1800 
frames  were  “read"  on  a  film  reader.  Film  reading  consists  of  setting  moveable 
cross-hairs  at  two  places  on  the  Image  of  the  model's  edge.  Hie  arc-tangent  pro¬ 
vides  the  angle  of  this  edge;  removal  of  the  cone  half-angle  then  gives  the  angle 
of  the  centerline  of  the  model. 

It  will  be  recalled  In  equations  (15)  or  (19)  that  the  angles  of  attack  and 
sideslip  must  be  available  at  the  same  time,  t-t^.  Thus  it  is  necessary  to 
designate  one  angular  record,  say,  a-a(t),  as  the  primary  and  interpolate  the 
second  record,  g"g(t)  In  this  case,  to  obtain  both  records  at  die  same  values  of 
time.  The  record  of  o(t)  (primary)  and  g(t)  (Interpolated)  are  given  in  Figure  9. 
in  alternate  presentation  of  essentially  the  same  information  Is  provided  In 
Figure  10.  Bara  a  (vertical)  is  plotted  with  g  (horizontal)  to  give  the 
familiar  a-g  petal  graphs.  Hie  damped  character  of  the  motion  is  more  obvious  in 
Figure  10  than  It  perhaps  la  In  Figure  9* 
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ANGLE  OF  ATTACK  - 

ANGLE  OF  SIDESLIP - 

TIME:  10  TO  16  SECONDS 
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ANGLE  OF  ATTACK  _ 

ANGLE  OF  SIDESLIP - 

TIME:  16  TO  20  SECONDS 


•  TIME  :  0  SECOND* 

•  TIME:  2.18  SECONDS 


ANGLE  OF  SIDESLIP  (DEGREES) 


FIGURE  10  (a)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDESLIP  FOR  TIME  INTERVAL  0  TO  2. IS  SECONDS 
(HCMMPHSRICAL  NOSE) 


®  TIME:  6.29  SECONDS 
®  TIME:  8.39  SECONDS 


ANGLE  OP  SIDESLIP  (DEGREES) 


m  ANGLE  OP  ATTACK  VtlltUS  ANGLE  OP  SIDESLIP  POR  TIME  INTERVAL  6.29  TO  &39  SECONDS 
pMNPmlnlvAL  NOvCf 


#  TIME:  8.30  SECONDS 
$  TIME:  10.54  SECONDS 


ANGLE  OF  SIDESLIP  (DEGREES) 


FIGURE  10  (•)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDESLIP  FOR  TIME  INTERVAL  8.30  T0 10.84  SECONDS 
(HEMI4PHERICAL  NOSE) 


%  TIME:  16.80  SECONDS 
®  TIME:  19.00  SECONDS 


ANGLE  OF  SIDE8LIR  (DEGREES) 


FIGURE  10  ill)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDESLIP  FOR  TIME  INTERVAL  16.80  TO  19.00  SECONDS 
(HEMISPHERICAL  NOSE) 


©  TIME:  14.72  SECONDS 
®  TIME:  16.80  SECONDS 


ANGLE  OF  SIDESLIP  (DEGREES) 


FIGURE  10  (I)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDESLIP  FOR  TIME  INTERVAL  14.72  TO  1S.80  SECONDS 
(HEMISPHERICAL  NOSE) 
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The  records  o(t^)  and  8(t^)  were  read  into  the  conputer  In  line  350  (see 
Table  2.  The  program  was  iterated  50  times.  The  last  20  iterative  steps  are 
given  in  Table  5. 

A  few  remarks  may  be  made  concerning  the  final  values  from  the  data  reduction 
program.  From  the  discussion  in  Section  3,  it  will  be  recalled  that  the  constants 
Alt  ®1*  A2*  ®2  are  of  no  value  *n  obtaining  the  aerodynamic  stability  derivatives. 
From  equation  (14)  it  is  possible  to  calculate  trim  angle.  Applying  this  equation 
to  the  values  of  and  of  Table  5  it  is  possible  to  calculate  the  trim  angle, 
esj,  as  1.12  degrees.  Initially  this  is  surprising  since  the  NAN  wind  tunnel 
model  does  not  have  any  apparent  configurational  asymmetries.  The  apparent  trim 
angle  is  nothing  more  than  an  error  in  the  background  fiducial  marking  used  for 
frame  alignment. 

The  frequency  and  damping  constants,  VI,  W2 ,  LI,  L2  are  of  value  in  computing 
the  three  aerodynamic  derivatives  and  the  spin  rata.  Inserting  these  values  as 
appropriate  in  the  program  of  Table  3  gives  the  derivatives  presented  in 
Table  4. 

The  static  pitching  moment  derivative  Cfflfl*‘CMA  of  -.328  is  the  value  of  this 
derivative  applicable  to  a  non-ablating,  spinning  conical  model  (7-degree  half 
angle,  0.223  nose  bluntness)  undergoing  bi-planar  angular  motion.  The  location  of 
the  statie  center  of  presaure  aft  of  the  nose  (fractions  of  body  length,  L)  is 
given  by  the  following: 

r  -  tt*  r  +  °*595  08) 

Is  L 

"o 

With  the  body  length  of  19.79  and  body  diameter  of  6.054  inches,  C  as  -.328  and 

“a 

Cg  of  1.35  the  center  of  pressure  is  0.669  body  lengths  aft  of  the  nose  and  0.0743 
o 

body  lengths  aft  of  the  c.g.  There  is  no  published  data  applicable  to  this  configu¬ 
ration  at  a  Mach  number  of  7.95.  However,  reproduced  from  an  unclassified  source  is 
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TABLE  5 

DATA  REDUCTION  FOR  DATA  POINTS  1  TO  900  FOR  SPHERE  CONE  MODEL. 


STANS. IEV.-E 

PREC-FRE0 

NU-FftEO 

PREC-0ANP 

NU-DANP 

.294289 

10.3401 

-14.2185 

-2. 4533 IE-2 

-2.06467E-2 

.294178 

10.3401 

-14.2175 

-2.42024E-2 

-2.10517E-2 

.294124 

10.3402 

-14.2147 

-2.59474E-2 

-2.12998E-2 

.294097 

10.3402 

-14.2142 

-2.S80O2E-2 

-2.14477E-2 

.294084 

10.3402 

-14.2158 

-2.54809E-2 

-2.15323E-2 

.294078 

10.3401 

-14.2155 

-2.55957E-2 

-2.I5779E-2 

.294074 

10.3401 

-14.2153 

-2.55348E-2 

-2.15999E-2 

.294073 

10.3401 

-14.2152 

-2.54913E-2 

-2.14081E-2 

.294072 

10.3401 

-14.2151 

-.02546 

-2.16085E-2 

.294072 

10.3401 

-14.215 

-2.54377E-2 

-.021605 

.294071 

10.3401 

-14.215 

-2.S4216E-2 

-2.15997E-2 

.294071 

10.3401 

-14.2149 

-2.54101E-2 

-2. 15939E-2 

.294071 

10.3401 

-14.2149 

-2.54018E-2 

-2.15883E-2 

.294071 

10.3401 

-14.2149 

-2.53958E-2 

-2.15833E-2 

.294071 

10.3401 

-14.2149 

-2.53915E-2 

-2.15789E-2 

.294071 

10.3401 

-14.2149 

-2.53884E-2 

-2.15752E-2 

.294071 

10.3401 

-14.2149 

-2.53862E-2 

-2.I5722E-2 

.294071 

10.3401 

-14.2149 

-2.53846E-2 

-2.15497E-2 

.294071 

10.3401 

-14.2148 

-2.53834E-2 

-2.15477E-2 

.294071 

10.3401 

-14.2148 

-2.53826E-2 

-2. 1544 IE-2 

Al 

11 

A2 

12 

A3 

.328411 

1 .74525 

.441837 

1.12199 

.017397 

B3 

lit 

112 

LI 

L2 

-9.10809E-3 

(0.3401 

-14.2148 

-2.53826E-2 

-2.15441E-2 

TABLES 

NON-DIMENSIONAL  AERODYNAMIC  COEFFICIENTS  AND  SPIN  RATE  FOR  DATA  POINTS  1  TO  900. 

(SPHERE-CONE  MODEL) 


STATIC 

CM 


IYMMIC 

CHS 


M9HUS 

CAPA 


SPIN 

P(RAS/8EC) 


C«P« 

xieler-LENSTH) 


* 


-.328249 


•1.4442 


.140737 


29.9489 


7.437071-2 
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data  for  this  configuration  at  Mach  10,  These  data  have  been  reproduced  in 
Figure  11 . 

It  will  be  noted  that  the  furthest  aft  that  the  center  of  pressure  reaches 
is  about  0.67  body  lengths.  As  just  noted ,  the  center  of  pressure  measured  in 
these  dynamic  tests  is  at  0.669.  It  is  not  possible  to  compare  the  static 
(constrained)  test  results  with  the  dynamic  test  results  since  the  angle  of 
attack  was  continuously  changing  during  the  dynamic  tests.  Nevertheless  the 
dynamic  testing  took  place  over  an  angle  of  attack  range  of  about  3  degrees. 

Taking  an  r.m.s.  value  of  angle  of  attack  at  something  like  2  degrees,  it  may 
be  seen  that  correlation  of  the  center  of  pressure  between  the  two  testing 
techniques  is  quite  satisfactory. 

4.2  ABLATING  NOSE  (AN)  WIND  TUNNEL  TESTS.  The  Three-Degree-of-Freedom  Gas 
Bearing  may  be  used  to  study  the  dynamic  motion  of  a  model  undergoing  surface 
ablation.  Ablation  was  confined  to  the  hemispherical  nose  by  replacing  the  steel 
nose  by  a  camphor  hemispherical  nose  molded  over  a  steel  bi-conic  substructure. 
This  substructure  is  shown  in  Figure  12.  The  particular  shape  of  this  sub¬ 
structure  was  derived  from  earlier  experimentation,  reference  (6). 

In  these  wind  tunnel  tests  the  low  temperature-ablating  nose  was  attached  to 
the  model.  The  testing  procedure  followed  along  the  lines  described  in  sections 
4.0  and  4.4.  Photographs  were  made  of  the  ablating  nose  throughout  the  20  seconds 
duration  of  the  erosion.  Eight  sequences  in  the  ablation  process  are  shown  in 


^Ragsdale,  W.C.,  at  al,  "Mach  10  IAP  Static  Force  Test  Program  in  the  NSWC/ 
WOL  Hypersonic  Tunnel  (WTR  1296),"  NSWC/W0L  MP  77-33,  September  1977. 

**  Jobe,  M.D.,  at  al,  "Bl-Modal  Flow  Field  Feasibility  Demonstration  in  the 
N8WC/H0  Hypersonic  Tunnel  (WTR  1319),"  NSWC  MP  79-217,  May  1979. 


M 

ay 


FIGURE  12  ABLATING  NOSE  SUBSTRUCTURE 
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Figure  13.  Each  photograph  depicts  phenomenon  about  2.5  seconds  apart.  Although 
not  evident  in  these  photographs,  the  model  is  under**'  simultaneous  pitch,  yaw 
and  spin.  It  may  be  seen  that  the  ablating  shape  approaches  the  bi-conic  sub¬ 
structure  illustrated  in  Figure  12. 

The  significant  difference  between  the  sphere-cone  model  and  the  ablating 
model  is  the  continuously  changing  nose  shape  associated  with  the  ablating  model. 
Changing  configuration  means  that  the  various  aerodynamic  derivatives  are  also 
continuously  changing  during  the  ablation  process. 

As  in  the  case  of  the  previously  discussed  sphere-cone  test,  primary  data 
consists  of  the  triplets  (a^,  61*  t^)  where  and  are  the  angles  of  attack 
and  sideslip  measured  at  time  t^.  Probably  the  most  straight  forward  way  of 
presenting  (a^,  8^,  t^)  data  points  is  as  a  time  history.  This  has  been  done  in 
Figure  14.  As  pointed  out  earlier,  the  two,  orthogonal  cameras  are  not  synchro¬ 
nized.  Since  there  are  some  differences  in  framing  rate  between  the  two  cameras, 
the  records  ot^  and  6^  are  not  at  the  same  time,  t^.  The  data  reduction  program 
does  require  the  two  angle  pairs  be  at  the  same  value  of  time.  Thus,  it  is 
necessary  to  Interpolate  one  angle  record  to  bring  it  into  coincidence  with  the 
other.  In  the  angular  history  presented  in  Figures  14,  the  sideslip  record, 
is  interpolated  to  coincide  with  the  angle  of  attack  record,  a^. 

Figures  15  are  an  alternate  form  of  representing  a^,  8^.  These  "petal" 
plots  are  given  at  about  2.0  second  intervals.  The  start  of  each  record  is 
Indicated  by  the  symbol Q  and  the  end  by  the  symbol^.  A  quick  perusal  of 
Figures  15  Indicates  slight  damping  until  about  10  seconds  (Figure  15f). 

The  succeeding  figures  seem  to  indicate  a  slight  increase  of  the  angular 
displacement  magnitude.  In  the  reduced  data  one  would  expect  to  find  changes 
in  sign  of  the  damping-ln-pltch  derivative,  C  . 
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FIGURE  13  NOSE  ABLATION  SEQUENCE 
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FIGURE  13  (CONTINUED) 
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In  the  first  reduction  of  these  data  the  whole  record  of  about  750  data  points 

was  divided  Into  sub- intervals  of  about  50  data  points  each.  Each  data  point  is 

characterized  by  the  triplets,  (a^,  8^»  t^) ,  i.e.  angle  of  attack,  angle  of  side 

slip  and  time.  A  data  record  of  50  data  points  has  a  time  duration  of  about  one 

second  as  has  been  pointed  out.  Thus,  each  petal  figure  consists  of  about  100 

data  points.  An  examination  of  the  aerodynamic  derivatives  reduced  from  these 

successive,  non-overlapping  sub-records  did  indicate  a  fairly  consistent  monotonic 

reduction  with  time  of  the  static  moment  derivative,  C  .  If  a  constant  value  is 

ma 

assumed  for  the  normal  force  derivative,  ,  the  preceding  observation  is 

a 

equivalent  to  recognizing  that  the  center  of  pressure  moves  steadily  forward  with 

time.  The  equivalent  statement  is  that  the  center  of  pressure  moves  forward  as 

the  model  ablates  from  a  sphere  cone  to  the  symmetrical  substructure  illustrated 

in  Figure  12.  As  will  be  seen,  the  center  of  pressure  moves  from  a  point  0.073 

body  lengths  (or  so)  aft  of  the  c.g.  to  a  point  0.054  body  lengths  aft  of  the  c.g. 

As  far  as  dynamic  measurements  are  concerned  it  is  felt  that  less  than  three 

significant  figure  accuracy  in  measurements  from  the  film  record  combined  with 

low  dynamic  pressure  compromised  reliable  dynamic  data.  In  addition  there  is  an 

indication  from  the  petal  plots  of  Figures  15  that  there  are  regions  of  both 

damped  and  undamped  motion.  From  the  reduction  of  the  50-polnt  records  the 

damping  in  pitch  derivative,  C  ,  varied  considerably  in  both  magnitude  and  sign. 

m 

q 

Sign  changes  (from  negative  to  positive)  can  be  used  to  qualitatively  indicate 
periods  of  damped  and  undamped  motion. 

In  order  to  examine  the  dynamic  measurements  in  more  detail  it  was  decided 
to  apply  the  data  reduction  method  to  overlapping  records,  still  50  data  points 
in  length.  This  examination  using  overlapping  records  was  confined  to  the 
beginning  of  the  ablation  record.  Thus,  records  of  data  points  1-50,  5-55, 

10-60,  etc.,  were  reduced.  The  results  of  this  data  reduction,  both  the 
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ANGLE  OF  ATTACK - 

ANGLE  OF  SIDESLIP - 

TIME:  5  TO  10  SECONDS 


•  TIME:  0  SECONDS 
$  TIME:  2.10  SECONDS 


m 

m 

«o 

uS. 


ANGLE  OP  SIDESLIP  (DEGREES) 


FIGURE  IB  <•)  ANGLE  OF  ATTACK  VERSU8  ANGLE  OF  SIDE  SLIP  FOR  TIME  INTERVAL  0  TO  £10 SECONDS 
(ABLATING  NOSE) 


•  TIME:  2.10  SECONDS 
0  TIME:  4.22  SECONDS 


ANGLE  Of  SIDESLIP  (DEGREES) 


FIGURE  IS  (b)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDE  SLIP  FOR  TIME  INTERVAL  110  TO  4.22  SECONDS 
(ABLATING  NOSE) 


ANGLE  OF  ATTACK  (DEGREES) 
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FIGURE  15  (c)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDE  SLIP  FOR  TIME  INTERVAL  4.22  TO  6.28  SECONDS 
(ABLATING  NOSE) 
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€>  TIME:  8.28  SECONDS 
$  TIME:  8.37  SECONDS 


ANQLE  OF  SIDESLIP  (DEGREES) 


FIGURE  18  (d)  ANQLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDE  SLIP  FOR  TIME  INTERVAL  6.28  TO  8.37  SECONDS 
(ABLATING  NOSE) 


W  TIME:  8.37  SECONDS 
®  TIME:  10.47  SECONDS 


Vk 


T1NE<  701)  >  14.0500 
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RUN  10 


n 


9  TIME:  14.08 SECONDS 
®  T»»E:  10.78  SECONDS 


ANQLfi  OF  SIDESLIP  (DIGRESS) 


Wa4,,l,  11  V8WUi ANOL0I OF stl DC SUE FOR TIME INTERVAL  14JB T0 10.78 S6CONDS 


r  <« *-H' *?v *••>;-< / 


ANGLE  OF  ATTACK  (DEGREES) 
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FIGURE  15  (i)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDE  SLIP  FOR  TIME  INTERVAL  16.75  TO  18.85  SECONDS 
(ABLATING  NOSE) 
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overlapping  and  non-overlapping  records,  la  given  in  Table  7.  The  number  of 
the  first  and  last  data  point  in  each  record  is  given  in  the  first  column.  The 
corresponding  values  of  tine  are  given  for  the  first  and  last  point  in  the  next 
column,  followed  by  the  average  time  for  the  record.  The  remaining  columns  give 

the  static  pitching  moment  derivative,  C  ,  the  pitch-damping  derivative,  C  ,  and 

m 

a  q 

the  center  of  pressure  aft  of  the  c.g.  in  body  lengths,  X.  It  is  felt  that  the 

spin  rate  was  too  low  for  the  Magnus  derivative  to  be  assigned  any  reliability; 

hence,  this  derivative  is  not  presented.  The  center  of  pressure  is  calculated 

from  the  pitching  moment  derivative,  C  ,  and  ah  estimated  value  of  the  normal 

ao 

force  derivative,  as  was  done  in  equation  (38).  The  difference  between  the  center 
of  pressure  defined  in  equation  (38)  and  in  aquation  (39),  below,  is  that  the 
reference  point  now  is  taken  as  the  center  of  gravity,  0.595  body  lengths  from 
the  nose. 


X  - 


\  (d/1) 


(39) 


A  value  of  C„  of  1.35  /  radian  was  taken  from  reliable  sources.  The  conversion 

a  _ 

(d/1)  is  used  to  express  X  In  terms  of  body  length  where  C  is  based  upon  the 

*o 

body's  maximum  diameter  as  a  reference. 

It  will  be  notad  in  Table  7.  that  the  center  of  pressure  at  the  beginning 
of  ablation  is  about  0.0724  body  lengths  aft  of  the  c.g.  This  value  might  be 
compared  with  the  value  of  0.0724  given  in  Table  6  for  the  sphere-cone  test. 

The  theoretical  value  for  the  non-ablating  sphere  cone  is  something  like  0.053 
body  lengths,  although  the  data  in  Table  7  is  for  an  ablating  model.  An 
examination  of  the  data  in  Table  7  shows  that  during  ablation  there  is  a  steady 
forward  motion  of  the  center  of  pressure.  At  the  beginning  of  ablation  the  center 
of  pressure  is  about  0.074  body  lengths  aft  of  the  c.g.  After  three  seconds, 
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TABLE  7  SUMMARY  OF  REDUCED  DATA  FROM  ABLATING 
MODEL  WITH  SYMMETRICAL  NOSE  SUBSTRUCTURE 


DATA  POINT 
NUMBERS 

TIME 

(SECONDS) 

AV.-TIME 

(SECONDS) 

C 

"a 

C 

"a 

X 

a  D 

1 

0 

50 

1.04 

0.52 

-.319 

-7.5 

0.0724 

5 

0.09 

55 

1.15 

0.62 

-.321 

-15.6 

0.0728 

10 

0.19 

60 

1.25 

0.72 

-.331 

-19.9 

0.0750 

15 

0.29 

65 

1.34 

0.81 

-.336 

“7.5 

0.0760 

20 

0.40 

70 

1.46 

0.93 

-.325 

+6.3 

0.0736 

25 

0.51 

75 

1.56 

1.035 

-.308 

+7.1 

0.0699 

50 

1.04 

100 

2.08 

1.560 

-.313 

+23.4 

0.0709 

100 

2.08 

150 

3.14 

2.61 

-.323 

-5.4 

0.0732 

150 

3.14 

200 

4.18 

3.66 

-.316 

-3.6 

0.0716 

200 

4.18 

250 

5.22 

4.70 

-.315 

11.0 

0.0714 

250 

5.22 

300 

6.26 

5.74 

-.314 

8.0 

0.0712 

300 

6.26 

350 

7.30 

6.78 

-.313 

-10.5 

0.0708 

350 

7.30 

400 

8.34 

7.82 

-.314 

-7.6 

0.0711 

400 

8.34 

450 

9.41 

8.87 

-.314 

+12.5 

0.0713 

450 

9.41 

500 

10.44 

9.92 

-.315 

+0.3 

0.0714 
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500 

10.44 

550 

11.49 

10.96 

-.315 

-11.5 

0.0714 

550 

11.49 

600 

12.50 

11.99 

-.315 

+4.0 

0.0713 

600 

12.50 

650 

13.59 

13.04 

i 

• 

N> 

00 

+7.1 

0.0650 

650 

13.59 

700 

14.63 

14.11 

-.241 

+1.7 

0.0545 

-»  .  -»■  --J*  ,(  y*  «  _  f  *  ■»  «  *-A  **»  •■  *  M. 
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the  center  of  pressure  Is  about  0.070  body  lengths  aft.  As  ablation  continues 
end  as  the  record  nears  the  end  the  center  of  pressure  moves  to  about  0.054  body 
lengths.  The  center  of  pressure  measured  in  a  wind  tunnel  on  a  constrained 
(etatic)  model  is  close  to  the  theoretical  value  of  0.053.  It  should  be  remembered 
that  the  model  in  the  dynamic  test  is  undergoing  motion  in  three  degrees  of  freedom 
and  that  the  windward  meridian  is  constantly  changing  in  body  coordinates. 

As  pointed  out  earlier,  the  damping  in  pitch  derivative,  C  ,  may  be  given 

mq 

only  a  qualitative  interpretation.  It  will  be  noted  at  the  beginning  of  ablation 
that  the  angular  record  indicates  dynamic  stability.  At  about  1  second  after  the 
start  of  the  record  the  motion  has  become  dynamically  unstable  until  something 
like  2  seconds  has  elapsed.  After  2  seconds,  the  motion  becomes  stable  again. 

This  alternate  dynamic  stability  And  instability  continues  throughout  the  record. 

Of  minor  interest  perhaps  is  the  fact  that  the  plane  of  oscillation  is 
rotating.  (See  Figures  15).  If  $  is  taken  as  the  rotation  rate  of  the  plane 
of  oscillation,  then  it  may  be  shown  that  there  is  a  relationship  between  ^  and 
the  nutatlonal  and  precessional  frequencies,  u>^  and  b^,  as: 


♦  ■  +  w2  (40) 

From  equation  (11c)  it  is  immediately  obvious  that : 


Therefore,  the  rate  of  rotation  of  the  plane  of  oscillation  is  a  simple  function 
of  the  roll  rate,  p. 

Of  considerably  more  interest  is  the  attitude  of  the  nose  toward  the  velocity 
vector  which  is  assumed  to  be  at  the  origin  of  the  a-0  petal  graphs.  It  will  be 
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ANGLE  OF  SIDESLIP  (DEGREES) 


FIGURE  IE  RELATIVE  POSITION  OF  INITIAL  WINDWARD  MERIDIAN  (  •  )  AND  INSTANTANEOUS 
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presumed  that  the  meridianal  line  on  the  body  (or  nose)  that  faces  the  origin  or 
velocity  vector  will  be  the  meridian  of  greatest  heat  transfer.  At  the  nose  this 
is  the  meridian  of  greatest  ablation.  Figure  16  is  a  rough  attempt  to  track 
the  windward  meridian  and  a  fixed  point  on  the  nose  as  the  model  undergoes  about 
one  oscillation.  The  fixed  point  on  the  nose  that  is  being  tracked  is  the 
original  stagnation  point  at  time  t  =*  0.00  seconds.  It  will  be  noted  that  the 
model  does  not  present  the  same  meridian  to  the  velocity  vector  as  the  model 
oscillates.  One  might  expect  fairly  uniform  ablation  of  the  model,  and  hence 
no  large  trim  angles,  should  develop  throughout  the  ablation  phase. 

It  would  seem  that  since  the  ablation  must  perforce  be  confined  to  the 
nose,  the  forces  (and  their  moments)  resulting  from  the  ablation  are  an  entirely 
localized  phenomenon.  This  may  not  be  entirely  accurate.  It  is  true  that 
there  are  local  forces  generated  by  the  local  ablative  mass  addition  to  the 
boundary  layer.  In  addition  however,  we  might  expect  the  generation  of  forces 
considerably  downstream  due  to  cross-flow  effects  because  the  boundary  layer 
has  been  thickened  by  the  upstream  ablation,  i.e.,  there  may  be  a  transport 
effect  associated  with  ablation.  In  addition  there  are  the  obvious  effects  due 
to  changing  nody  geometry  resulting  from  the  ablation.  In  the  tests  under 
discussion,  with  ablation  confined  to  the  nose,  body  shape  changes  are  limited 
to  the  Increase  in  nose  bluntness  and  the  development  of  a  bi-conic  shape  (see 
Figures  12  and  13).  Generally  an  Increase  in  nose  bluntness  causes  a 
decrease  in  the  damping  of  angular  motion.  The  effects  caused  by  changes  in 
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injection  into  the  boundary  layer,  either  locally  or  by  transport  through  the 
boundary  layer  from  some  other  part  of  the  body."  ^ 

The  question  still  remains  as  to  whether  or  not  ablation  situated  at  the 
model's  nose  causes  dynamic  instability.  The  test  results  reported  herein  in¬ 
dicate  that  nose  ablation  does  cause  dynamic  instability.  These  results  agree 
with  the  conclusions  of  Ericsson  and  Reding  (reference  (5)): 

"A  slender  conical  reentry  body  may  experience  increased  static,  stability 
but  severely  decreased  dynamic  stability  at  low  angles  of  attack  and  low 
oscillation  amplitudes  when  ablation  is  concentrated  to  the  blunted  nose 
and  cone  trailing  edge." 

Ericsson  and  Reding' s  paper  (reference  (5))  is  an  analytic  exercise  and  is  not 
substantiated  by  any  experimental  results.  The  authors  admit  in  the  conclusion 
to  their  paper:  "...  the  analysis  presented  herein  is  highly  speculative  ...". 
Nevertheless  to  some  extent  the  conclusions  of  this  paper  as  indicated  by  the 
above  quotation  seem  to  be  substantiated  by  the  results  of  the  ablating  nose 
test  reported  herein. 

In  1965  a  series  of  experiments  were  carried  out  to  examine  the  dynamic 
behavior  of  an  ablating  model  undergoing  single  degree  of  fre  iom  angular  (pitch) 
motion.  These  experimental  results  were  reported  in  what  may  be  regarded  as  the 
seminal  paper  in  the  experimental  investigations  of  the  effects  of  surface 

g 

ablation  on  the  dynamics  of  a  reentry  body.  While  this  paper  by  Grimes  is 
based  upon  wind  tunnel  dynamic  tests,  it  should  be  appreciated  that  there  are 

^Ericsson,  L.E.,  et  al,  "Ablation  Effects  on  Vehicle  Dynamics,"  Journal  of 
Spacecraft ,  Vol.  3,  No.  10,  October  1966. 

8Grimes,  J.H.,  "Influence  of  Ablation  on  the  Dynamics  of  Slender  Reentry 
Configuration,"  Journal  of  Spacecraft,  Vol  2,  Jan/Feb  1965. 
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considerable  differences  between  the  model  support  system  used  in  Grimes 
experiments  and  the  three-degree-of-freedom  gas  bearing.  The  model  in  the  tests 
reported  in  reference  (8)  cannot  spin  and  its  angular  motion  is  restricted 
entirely  to  a  plane.  By  way  of  contrast  the  model  in  the  three-degree-of-freedom 
gas  bearing  is  free,  within  an  upper  bound  restriction  on  total  angle  of  attack, 
to  undergo  simultaneous  and  mechanically-uncoupled  yaw,  pitch  and  roll  motion. 

In  other  words  a  comparison  is  being  made  between  experimental  results  from 
single  degree  of  freedom  and  three  degree  of  freedom  model  supports.  Clearly 
there  must  be  considerable  differences  between  the  ablation  physics  in  both 
cases.  In  the  three-degree-of-freedom  support  the  presence  of  roll  rate  allows 
the  region  of  ablation  to  alter  its  orientation  with  respect  to  the  velocity 
vector  independent  of  the  angle  of  attack.  Consequently  neither  Grimes  nor 
Ericsson's  paper  for  that  matter  can  be  entirely  supportive  of  the  gas  bearing 
test  presented  in  this  report.  Indeed  there  seems  to  be  some  conflict  between 
Ericsson  and  Grimes  as  to  the  effect  of  ablation  on  dynamic  stability.  The 
following  quotation  is  taken  from  Grimes'  paper  which  conflicts  with  the 
experimental  results  of  this  report  and  the  analytic  results  of  Ericsson  as 
Indicated  by  the  previously  given  quotation. 

"Coating  the  entire  model  ...  forward  of  the  center  of  gravity,  produced 
a  dynamically  stable  configuration  with  Cmq  values  approaching  those  for 
flight  test.  Coating  the  rearward  sections  of  a  conic  or  cone-cylinder 
flair  produces  a  dynamically  unstable  model." 

In  the  above  quotation  the  word  "coating"  refers  to  the  application  to  the  model 
of  a  low  temperature  ablator  such  as  paradlchlorobenzine  CgH^Clj*  Locating 
the  "coating"  is  equivalent  to  setting  the  location  of  the  model's  surface 


ablation 


In  brief  the  results  of  the  three-degree-of-freedom  support  indicate  that 
while  the  model  is  initially  dynamically  stable,  the  initiation  of  ablation 
causes  dynamic  instability  which  continues  until  the  ablation  ceases  and  the 
metalic  bi-conic  nose  is  revealed. 

The  following  is  a  partial  listing  of  the  conclusions  that  might  be  drawn 
from  these  ablation  tests: 

(a)  The  static  center  of  pressure  is  at  about  0.073  body  lengths  aft  of  the 
center  of  gravity  (0.595  body  lengths  aft  of  the  nose)  at  the  beginning  of 
ablation.  This  value  is  essentially  the  same  as  that  measured  for  the  steel 
nose  (non-ablating)  sphere  cone  model. 

(b)  The  static  center  of  pressure  moves  forward  during  ablation  to  a  value 
of  about  0.054  body  lengths. 

(c)  Agreement  seems  to  be  acceptable  with  results  from  supported  testing 

where  the  center  of  pressure  is  about  0.053  body  lengths.  It  must  be  remembered 

that  dynamic . testing  about  a  single  axis  cannot  provide  the  normal  force 

derivative,  ,  and  hence  a  final  evaluation  of  the  center  of  pressure  must  rely 
a 

on  an  external  source  of  C,t  . 

Na 

(d)  Dynamic  measurements  indicate  regions  of  dynamic  stability  and 
Instability.  No  convincing  argument  is  available  at  this  time  to  entirely 
explain  this  observation. 

4.3  POST-ABLATION  WIND  TUNNEL  TEST  -  ABLATING  NOSE  SUBSTRUCTURE 
(SYMMETRICAL) .  In  the  previous  section  (Section  4.2))  the  ablating  nose  tests 
were  discussed.  It  should  be  recalled  and  emphasized  that  ablation  was  con¬ 
fined  to  the  nose  by  replacing  the  steel  nose  by  a  low  temperature  ablator. 

The  inital  shape  of  this  ablating  nose  was  a  hemi-sphere  of  0.223  bluntness 
(see  Figure  7).  The  substructure  of  this  nose  is  the  contour  illustrated 
in  Figure  12.  When  the  ablation  process  is  completed  this  contour  is 
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revealed.  The  "post-ablation"  shape  was  developed  from  photographic  records 
of  preserved  contours  from  earlier  static  ablation  tests.  It  may  be  seen  in 
the  ablation  sequences  of  Figures  13  that  the  contour  does  ablate  to  something 
like  this  shape,  so  that  Figure  12  depicts  the  asymtotic  contour  of  the  ablation 
process.  It  is  of  course  necessary  to  preserve  this  ablated  shape  as  a  steel 
substructure  for  post-ablation  dynamic  testing  -  the  low  temperature  ablator 
would  continue  to  change  shape  under  ambient  conditions. 

It  was  assumed  that  60  seconds  after  initiation  of  the  ablation  run,  that 
ablation  was  complete  and  the  steel  substructure  was  completely  revealed.  The 
kicking  jet  was  then  Initiated  to  supply  angular  displacement  to  the  model. 

This  would  then  begin  a  new  run  designated  as  the  post  ablation  run  with 
symmetrical  steel  sub-structure.  Since  it  was  not  possible  to  be  completely 
certain  that  ablation  has  removed  all  of  the  low  temperature  ablator,  it  was 
decided  toreduce  the  750  data  point-long  run  in  overlapping  records  of  150 
data  points  each.  The  result  of  this  exercise  is  presented  in  Table  8.  It 
will  be  noted  that  the  center  of  pressure  is  nearly  exactly  that  measured  at 
the  end  of  the  ablating  run  (see  Table  7).  The  a-8  history  of  this  post¬ 
ablation  run  is  given  in  Figures  17;  the  a-8  petal  plots  are  presented  in 
Figures  18.  It  will  be  noted  that  these  petal  plots  for  the  post-ablation  run 
are  nearly  circular  as  compared  with  the  more  leaf-like  shape  observed  during 
the  ablation  run  (see  Figures  15)  as  the  sphere-cone  runs  (see  Figures  10). 

It  is  possible  to  show  that  the  character  of  the  motion  is  merely  dependent 
upon  initial  conditions  in  a  and  8.  First  by  way  of  definition. 
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SUMMARY  OF  REDUCED  DATA  FROM  POST-ABLATION 
WITH  SYMMETRICAL  NOSE  SUBSTRUCTURE 


DATA  POINT  TIME  AV.-TIME  C  C 

NUMBERS  (SECONDS) _ (SECONDS) _ “a _ j 
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ANGLE  OF  ATTACK - 

ANGLE  OF  SIDESLIP - 

TIME:  10  TO  15  SECONDS 


«>  TIME:  £09  SECONDS 
$  TIME:  4.17  SECONDS 


ANGLE  OF  SIDESLIP  (DEGREES) 


FIGURE  IS  (b)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDE  SLIP  FOR  TIME  INTERVAL  2.00  TO  4.17  SECONDS 
(ASUTINQ  NOSE  SUBSTRUCTURE  -  NO  TILT) 


ANGLE  OF  ATTACK 


ANGLE  OF  StDESUF  (DEGREES) 


V  n  ■  a 

t"0  c 
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£'“° 


H/t-°  ■  ®o 

In  equations  (6)  it  was  shown  that  the  motion  is  Initially  established  by 
the  8 lx  constants  A^,  B^,  A^,  B^,  A^,  By  One  of  the  constants,  A^,  was  given 
in  equations . (6)  as  a  generic  form.  Most  reentry  vehicles,  as  well  as  other 
statically  stable  vehicles,  have  the  characteristic  that  the  damping  exponents 
X^  and  X^  are  small  in  comparison  to  the  circular  frequencies,  and  as 


i 

i 


hK<  wi 


X2  «  o»2 


Neglecting  X ^  and  X2  where  they  appear  in  the  expressions  for  the  above  six 
constants  results  in  an  error  that  is  typically  about  2  percent.*  With  the 
assumption  noted  the  six  constants 


*  A  more  exact  estimate  of  this  error  is 


X100 


which  for  the  wind  tunnel  runs  under  consideration  are  about  2Z. 
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nay  be  written  as. 


k  .  *_°  *  “2  “o  +  (P  ~  M2>  A3  ~  K  *  m2  ao 

1  "  "  [  “l  -  “2  J  "  [  “l  - 


a  -  u>,  B 

,  .  _2 - 2 — o 

1  -  u>2 


S_  +  w,  a  +  (p  -  to.)  A  B  +»,  a 

O  1  O  1  3  iv  o  1  o 

A2  w2  “  W1  w2  -  wT 


“o  I  “l  go 
U2’U1 


v  r _ i  1 

ly  W2  (p  "  Ul)  +  p(ul  “  p) 


>  P  ^  tu. 


B3  “  0  ,  p  i  (Oj. 


where  the  simplification  at  the  far  right  assumes  no  body— fixed  asymmetry,  i.e.. 


A^  ~  <5  -  0 .  If  and  are  regarded  as  the  components  of  a  polar  vector, 
in  the  8  (horizontal)  and  a  (vertical)  direction,  respectively,  then. 


i2  .  _2 


*4.1  t-o  "  K10  "  tAl  +  B1 


(44a) 


A  similar  relationship  may  be  written  far  A2  and  B2  as. 


K2 1 1-0  K20  "  ^A 2  +  B2 


(44b) 


and  for  A„  and  as. 
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*3  W)  -  K30  '  VA3  "  B3  <4 

It  can  easily  be  shown  from  equations  (6)  and  (10)  the  rotates  with  a  rate 

radious/ second  in  the  same  direction  as  the  roll  rate,  p  and  that  rotates 

at  a  rate  in  a  direction  opposite  to  the  spin  rate  K^,  which  accounts  for  a 

body-fixed  asymmetry  rotates  at  the  spin  rate,  p. 

The  geometry  of  this  "tri-cydic  motion"  is  given  in  the  sketch  below. 


MODEL  VERTEX 


a 


ANGLE  OF  SIDEtUF  (DEGREES) 


FIGURE  19  (a)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDESLIF  FOR  TIME  INTERVAL  0  TO  IBS  SECONDS 
(ABLATING  NOSE  SUBSTRUCTURE  -  3  DEG  TILT) 


ANGLE  OF  SIDESLIP  (DEGREES) 


FIGURE  IS  (b)  ANGLE  OF  ATTACK  VERSUS  ANGLE  OF  SIDESLIP  FOR  TIME  INTERVAL  Z08  TO  4.16  SECONDS 
(ABLATING  NOSE  SUBSTRUCTURE  -  3  DEG  TILT) 
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In  examining  the  motion  of  the  model  in  Figures  18,  the  observed  circular  motion 

is  due  to  a  vanishing  of  the  Kg  arm.  Now,  it  will  be  recalled  that  the  air- jet 

acts  only  to  cause  an  angular  displacement  in  the  vertical  plane,  i.e.,  to  set 

a  (and  perhaps  a  ).  However,  since  the  post-ablation  run  was  initiated  with 
o  o 

some  residual  motion  from  the  preceding  ablating  run,  it  is  likely  that  6 

o 

and  8q  were  not  zero  at  initiation  of  the  kicking  jet. 

A  fortuitous  set  of  circumstances  could  well  result  in  the  minimisation  or 
elimination  of  the  Kg  arm  by  reducing  or  eliminating  Ag  and/or  Bg.  This  could 
happen  if  the  numerators  in  the  expressions  for  Ag  and  Bg  in  equations  (43)  are 
small.  For  example,  remembering  that  Wg  is  always  negative  (p,  spin  rate,  and 
nutation  rate  taken  as  positive),  if  is  negative  and  aQ  is  positive,  then 
if 

$o  "  '  “l  “o  C45) 

the  numerator  of  Ag  vanishes.  The  operation  of  the  kicking  jet  also  occurs 
when  both  and  are  small,  reducing  the  magnitude  of  Bg  (as  well  as  Ag). 

The  point  of  this  discussion  is  that  the  presence  of  circular  motion  does 
not  indicate  any  unusual  dynamic  effects.  A  comparison  of  the  center  of  pressure 
values  between  the  ablation  test  (Table  7)  and  the  post-ablation  test  (Table  8) 
indicates  coincidence  in  this  parameter. 

4.4  POST-ABLATION  WIND  TUNNEL  -  ABLATING  NOSE  SUBSTRUCTURE  (ASYMMETRICAL). 
The  test  discussed  in  this  section,  like  that  of  the  previous  section,  is  also 
a  test  of  the  model  with  the  ablating  nose  substructure  in  place  (see  Figure  12). 
However,  unlike  the  previous  test  the  axis  of  symmetry  of  this  nose  is  tilted 
3  degrees  with  respect  to  the  center  line  of  the  conical  after  body.  This  nose 
tilt  then  Introduces  an  asymmetry  into  the  configuration.  The  results  of  this 
test  are  given  below  in  Table  9.  The  a-8  petal  plats  are  given  in  Figure  19. 
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-.2575? 

-1.98461 
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-33.4564 
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It  will  be  noted  in  Table  9  that  the  center  of  pressure  is  comparable  to  that 
given  for  the  symmetrical  ablating  nose  substructure  (Table  8).  Figures  19 
do  not  indicate  the  circular  motion  observed  for  the  symmetrical  nose.  The 
asymmetric-nose  tests  were  performed  as  an  Initial  inspection  in  the  wind  tunnel 
and  did  not  follow  an  ablating  test.  Consequently  and  $Q  were  negligibly 
small  so  both  the  and  Kj  modes  were  excited. 

5.0  CONCLUSIONS 

These  tests  led  to  some  interesting  results.  It  was  shown  that  the  no- - 
ablating  sphere-cone  model  had  a  center-of-pressure  that  was  0.669  body  lengths 
aft  of  the  center  of  gravity.  Comparisons  were  made  favorably  with  earlier 
static  wind  tunnel  data.  The  ablating  nose  model  had  a  center  of  pressure 
almost  exactly  at  the  same  point. as  the  non-ablating  model  at  the  on  set  of 
ablation.  During  ablation  there  waa  a  steady  foward  motion  of  the  center  of 
pressure.  Dynamic  measurements  during  the  ablation  indicate  that  there  are 
regions  of  dynamic  stability  and  instability. 

Completion  of  this  test  program  has  been  more  laborious  than  anticipated 
due  to  the  failure  of  the  fiber  optic  readout  system  to  function  satisfactorily. 
Nevertheless,  it  has  been  demonstrated  that  the  three-degree-of-freedom  model 
support  can  be  used  to  make  pitching  moment  measurements  on  ablating  and  non- 
ablating  models  in  the  Naval  Surface  Weapons  Center's  Hypersonic  Wind  Tunnel. 
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